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TERRAIN  BACKGROUND  SIGNATURE  RESEARCH 


PART  I :  INTRODUCTION 

Background 

1.  The  detection  of  camouflaged  targets  that  are  located  in  a  terrain 
background  generally  involves  a  comparison  of  the  target  signature  and  the 
terrain  background  signature.  This  is  the  case  for  active  and  passive  sensors 
and  for  all  frequency  bands  in  which  the  sensor  operates- -microwave  (mw) , 
millimeter  wave  (mmw) ,  infrared  (IR),  and  visible.  The  magnitude  of  the 
target-background  signature  contrast  is  used  for  the  prediction  of  target 
detection  probabilities,  for  false  alarm  rates,  and  for  any  general  estimate 
of  target  visibility.  Modern  technology  has  provided  sophisticated  devices 
that  offer  increased  signature  discrimination  capabilities  that  enhance  target 
detection  performance.  These  include  radars,  forward-looking  IR  sensors 
(FLIRS),  laser  scanners,  and  target  locators. 

2.  Two  main  environmental  factors  determine  the  target  location  per¬ 
formance  of  electromagnetic  sensors- -terrain  and  weather.  The  terrain  back¬ 
ground  environmental  factors  include  soil  type,  soil  moisture,  grass,  bushes, 
forest  canopies,  agricultural  development,  and  urban  development.  Targets  can 
be  camouflaged  so  as  to  blend  into  the  terrain  background  by  exhibiting  a 
signature  similar  to  the  terrain  background  signature  and  having  a  small 
target -background  signature  contrast.  Also,  targets  can  be  partially  or  com¬ 
pletely  hidden  from  view,  as  in  the  case  of  terrain  with  sufficiently  large 
mesoroughness  due  to  hills  or  vegetation  such  as  tall  grass,  bushes,  and  trees 
(Weiss  1981,  1982).  Weather  conditions  often  render  sophisticated  electro- 
optical  sensor  systems  ineffective.  This  may  be  the  result  of  smog,  haze, 
clouds,  rain,  and  snow  that  obscure  a  target  from  a  sensor  (Kays  et  al.  1980). 
Environmental  effects  predictions 

3.  Radar,  IR,  and  laser  target  detectors  are  affected  by  environmental 
conditions  in  their  own  unique  ways.  Therefore,  for  a  given  set  of  environ¬ 
mental  factors,  some  electromagnetic  sensors  will  operate  more  efficiently 
than  others,  and  predicting  which  system  will  perform  in  an  optimum  way 
becomes  necessary.  The  basis  of  these  predictions  is  a  complete  knowledge  of 
the  absorption,  emission,  scattering,  and  reflectance  of  electromagnetic  waves 
from  various  target  and  terrain  background  surfaces,  and  the  transmission  and 
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attenuation  of  electromagnetic  waves  through  the  atmosphere.  Only  in  this  way 
can  the  electromagnetic  target/background  signature  contrast  seen  by  a  distant 
sensor  be  predicted  in  terms  of  prevailing  weather  and  terrain  conditions. 

4.  A  combined  experimental  and  theoretical  investigation  of  the  physi¬ 
cal  phenomena  involved  in  the  detection  of  targets  located  in  complex  terrain 
backgrounds  has  been  undertaken  by  the  US  Army  Corps  of  Engineers .  This 
investigation  required  the  development  of  analytical  models  that  simulate  the 
physical  effects  of  the  interaction  of  electromagnetic  waves  with  the  sensor, 
atmosphere,  target,  and  terrain  background.  Thus,  the  effects  of  weather  and 
terrain  on  the  performance  of  electromagnetic  sensors  can  be  predicted  for  the 
time  of  engagement  with  the  enemy,  and  the  optimum  type  of  electromagnetic 
sensor  can  be  determined  for  the  expected  weather  and  terrain  conditions. 

5.  The  interaction  of  electromagnetic  radiation  with  the  weather  struc¬ 
tures  of  the  atmosphere  and  with  the  various  terrain  background  textures  con¬ 
sists  of  complicated  physical  processes  whose  description  requires  an  under¬ 
standing  of  the  basic  physics  involved  in  these  processes.  Therefore, 
physics-based  analytical  models  are  required  to  simulate  the  matter- radiation 
interaction  processes  associated  with  the  components  of  target -detection  per¬ 
formance  models.  These  models  are  necessarily  approximations  to  the  real 
effects  that  occur  in  the  radiation- environment  system,  and  the  user  must  be 
aware  of  the  limited  validity  of  the  physics  models.  The  user  must  also  be 
aware  of  the  extreme  variability  of  the  terrain  background  and  weather  struc¬ 
tures,  as  well  as  the  extreme  variability  of  the  associated  physical  param¬ 
eters  required  by  the  models  to  make  predictions.  The  main  reason  for  the 
discrepancies  between  model  predictions  of  target  and  terrain  background  sig¬ 
natures  and  the  actual  measured  target  and  terrain  background  signatures  is 
the  lack  of  accurate  knowledge  of  electromagnetic  and  thermal  parameters  of 
targets,  terrain  backgrounds,  and  weather  for  the  time  when  signature  measure¬ 
ments  are  obtained. 

6.  Many  modem  weapon  systems  operate  via  a  fusion  of  several  types  of 
sensors  so  that  it  becomes  important  to  develop  models  which  predict  the 
multispectral  appearance  of  a  target  in  a  background  and  calculate  the  attenu¬ 
ation  of  the  various  modes  of  electromagnetic  energy  transfer  from  the  target 
to  the  sensor.  Acoustic  and  seismic  modes  of  energy  transfer  are  also  pos¬ 
sible;  however,  this  report  considers  only  the  transfer  of  energy  by  electro¬ 
magnetic  waves  such  as  in  the  case  of  radar  (mmw  and  mw) ,  thermal  FL1RS 
(0.7100  pm) ,  and  laser  scanners  and  visible  light  cameras  (0.4-0. 7  pm). 
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Weather  conditions  affect  the  operation  of  all  these  devices,  and  the 
analytical  models  that  are  used  for  target  signature  prediction  must  be  multi - 
spectral,  weather - dependent  computer  codes. 

7.  The  target  detection  capability  of  a  sensor  is  often  expressed  as  a 
target  detection  probability  that  is  the  product  of  several  component  prob¬ 
abilities:  for  the  target/background,  for  the  atmospheric  conditions,  and  for 
the  electromagnetic  sensor.  The  prevailing  weather  conditions  affect  both  the 
target/background  target  detection  probability  and  the  atmospheric  component 
of  the  target  detection  probability.  For  instance,  the  IR  target/background 
contrast  is  minimized  on  a  rainy  day,  and  the  radar  returns  from  wet  vegeta¬ 
tion  often  mimic  the  radar  returns  from  metal  targets .  An  empirical  determi - 
nation  of  all  possible  environmental  influences  on  the  success  of  a  weapon 
sensor  system  is  very  costly  and  time  consuming;  consequently,  it  is  advanta¬ 
geous  to  use  analytical  models  to  predict  sensor  performance  for  prevailing 
weather  and  terrain  conditions. 

8.  The  physical  characteristics  of  terrain  backgrounds  such  as  physi¬ 
cal  temperature,  emissivity,  absorptivity,  heat  conductivity,  dielectric  con¬ 
stant,  magnetic  permeability,  and  soil  moisture  content  vary  with  the  time  of 
day  and  season  because  of  the  diurnal  and  seasonal  variation  of  the  insolation 
and  the  effects  of  local  weather  variations.  Therefore,  the  thermal,  radar, 
and  laser  signatures  of  terrain  backgrounds  will  change  with  weather  and  sea¬ 
son  in  addition  to  having  a  diurnal  variation.  Target  signatures  will  also 
experience  these  changes,  and  therefore  a  comprehensive  target/background 
signature  model  should  be  able  to  predict  the  multispectral  contrast  between  a 
target  and  background  by  simulating  the  effects  of  season,  weather,  and  diur¬ 
nal  solar  flux  variation. 

9.  The  models  and  algorithms  described  in  this  report  are  designed  to 
be  applied  to  the  problem  of  locating  vehicles  and  high-value  targets.  Of 
high  value  are  large  critical  facilities  such  as  airfields;  petroleum,  oil, 
and  lubricants  (POL)  storage  tanks;  bridges;  power  plants;  and  rail  yards. 

Each  target  will  have  a  characteristic  multispectral  signature  whose  contrast 
with  the  surrounding  terrain  background  signatures  must  be  predicted  at  the 
time  of  engagement  with  the  enemy. 

10.  Much  work  has  been  done  in  the  area  of  target -background  signature 
contrast  prediction  both  at  the  US  Army  Engineer  Waterways  Experiment  Station 
(WES)  and  elsewhere.  This  report  summarizes  the  WES  contribution  to  multi¬ 
spectral  target -background  signature  simulation  with  application  to  thermal 
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IR,  laser,  and  radar  target  locators.  The  physical  structure  and  input  param¬ 
eters  of  the  signature  prediction  models  are  described  with  consideration 
given  to  the  immense  complexity  of  the  problem  of  forecasting  the  performance 
of  sophisticated  weapon  systems. 

Framework 

11.  In  addition  to  a  sensor  characteristics  model  and  a  general  systems 
integration  model  that  combines  the  basic  computer  programs,  two  basic  compo¬ 
nents  are  needed  for  an  electromagnetic  sensor  performance  prediction  model: 

(a)  a  targe t/background  model  for  laser,  radar,  and  laser  sensor  systems  and 

(b)  an  electromagnetic  and  thermodynamic  parameter  database  for  targets  and 
terrain  backgrounds  that  will  be  used  to  drive  the  terrain  background  signa¬ 
ture  models. 

12.  Often  a  target  and  a  terrain  background  are  sufficiently  simple  so 
as  to  be  described  by  a  simple  signature  contrast  value  for  the  purpose  of 
target  location.  For  this  case,  the  differences  in  visible  light  reflec¬ 
tivity,  infrared  radiance,  and  radar  cross  section  can  be  used  to  discriminate 
targets  from  backgrounds.  If,  however,  the  target  and  terrain  background  are 
complex  and  contain  large  variations  of  reflectivity,  radiance,  and  radar 
cross  section,  local  "hot  spots"  must  be  considered  on  the  target  and  in  the 
terrain  background,  and  the  simple  concept  of  signature  contrast  becomes  inad¬ 
equate  to  predict  target  detectability.  A  statistical  analysis  of  the  spatial 
variation  of  the  background  signatures  may  be  of  value,  but  clearly  for  this 
case  a  full  scene  simulation  and  an  imaging  sensor  should  be  used. 

13.  The  signature  calculations  that  would  enter  a  full  scene  simulation 
for  laser,  radar,  and  thermal  IR  sensors  are  performed  by  describing  the 
unique  physical  processes  that  occur  when  electromagnetic  waves  in  these  char¬ 
acteristic  wavebands  interact  with  or  are  transformed  by  a  target  and  terrain 
background  surface.  For  laser  light  interacting  with  a  terrain  surface,  the 
relevant  processes  are  surface  scattering  and  reflectance,  and  a  calculation 
of  the  bidirectional  reflectance  is  required  in  terms  of  frequency,  depression 
angle  of  the  source,  and  view  angle  of  the  sensor.  For  radar,  the  relevant 
processes  are  scattering  and  reflection,  and  the  radar  differential  cross 
section  must  be  calculated  in  terms  of  surface  geometry  and  the  electrical 
parameters  of  the  terrain  surface.  The  calculation  of  the  IR  signature  of  a 
terrain  background  requires  the  evaluation  of  an  energy  budget  at  the  terrain 
surface,  which  includes  all  of  the  relevant  modes  of  energy  transfer  at  the 
surface.  This  determines  the  physical  temperature  and  radiance  of  the  terrain 
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surface.  The  modes  of  energy  transfer  are  determined  by  the  local  weather  and 
the  thermal  properties  of  the  ground. 

Environmental  parameters 
for  thermal  signature  models 

14.  The  primary  driving  factors  that  determine  the  thermal  IR  signa¬ 
tures  of  terrain  backgrounds  are  local  weather  and  type  of  terrain  surface, 
such  as  soil,  vegetation,  asphalt  runway,  etc.  The  weather  affects  the  amount 
of  solar  radiation  incident  on  the  terrain  surface  and  therefore  can  indi¬ 
rectly  affect  the  surface  conditions  such  as  water  content.  Weather  directly 
affects  the  surface  temperature  by  regulating  the  incident  solar  radiation  and 
sky  radiation,  which  is  transformed  to  heat  at  the  terrain  surface  (Balick  et 
al.  1981).  An  atmospheric  model  is  needed  to  determine  the  amount  of  attenua¬ 
tion  of  the  background  radiance  due  to  absorption  and  scattering  by  the 
atmosphere . 

15.  The  thermal  models  require  their  own  characteristic  terrain  back¬ 
ground  and  weather  information.  The  required  terrain  and  geographic  data 
include  terrain  elevations,  soil  moisture  content,  vegetation  moisture  con¬ 
tent,  heat  conductivity  and  diffusivity,  shortwave  absorptivity  and  longwave 
emissivity  of  the  terrain  surface,  latitude  of  area  of  interest,  and  the 
zenith  and  azimuth  angles  of  the  sun  (Balick  et  al.  1981;  Balick,  Scoggins, 
and  Link  1981).  The  weather  data  include  air  temperature,  cloud  cover,  cloud 
type,  relative  humidity,  wind  speed,  IR  sky  radiation,  and  insolation.  The 
weather  and  terrain  background  conditions  must  be  known  for  the  precise  time 
of  engagement  with  an  enemy  target.  This  is  the  primary  drawback  of  using 
analytical  models  for  thermal  IR  signature  prediction  because  the  terrain 
background  parameters  are  known  imprecisely  if  at  all. 

16.  A  prerequisite  for  the  determination  of  thermal  IR  images  of  tar¬ 
gets  and  backgrounds  is  the  prediction  of  target  and  background  temperatures. 
Many  target  and  terrain  temperature  estimation  models  have  appeared  in  the 
literature  (Kahle  1977;  Deardorff  1978;  Link  1979;  Jacobs  1980;  Smith  et  al. 
1980;  Balick  et  al.  1981;  Balick,  Scoggins,  and  Link  1981;  Botkin  et  al.  1981; 
Baird  et  al.  1982;  Hechinger,  Raffy,  and  Becker  1982;  Solomon  1982;  Hodge 
1989) .  The  WES  has  for  several  years  conducted  research  on  the  camouflage  of 
fixed  industrial  and  military  installations.  This  work  included  a  research 
program  for  developing  a  realistic  target  and  surround  database  for  use  in  the 
design  and  evaluation  of  imaging  and  nonimaging  sensors  for  surveillance, 
target  acquisition,  and  terminal  homing.  A  preliminary  procedure  for 
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predicting  terrain  surface  temperatures  was  developed  during  this  research 
(Link  1979). 

17.  The  prediction  of  terrain  surface  temperatures  uses  an  energy  bud¬ 
get  condition  at  the  ground  surface  (Kahle  1977,  Balick  et  al.  1981).  The 
effects  of  insolation,  atmospheric  IR  emission,  heat  conduction  and  convection 
in  the  air,  thermal  IR  emission  from  the  ground  surface,  heat  conduction  in 
the  ground,  and  latent  heat  of  evaporation  and  condensation  are  used  to  form 
the  energy  budget.  A  one -dimensional  (1-D)  analysis  using  this  energy  budget 
produced  the  Terrain  Surface  Temperature  Model  (TSTM) ,  which  applies  to  non- 
vegetative  surfaces  such  as  soil,  rock,  asphalt,  and  concrete  (Balick  et  al. 

1981)  .  Extensions  to  vegetative  canopies  have  been  considered  (Deardorff 
1978;  Smith  et  al.  1980;  Balick,  Scoggins,  and  Link  1981).  An  initial  two- 
dimensional  (2-D)  generalization  of  these  heat  flow  concepts  has  also  been 
developed  for  nonvegetated  surfaces  (Solomon  1982).  Several  calculations  of 
the  temperature  of  vertical  walls  have  appeared  in  the  literature  (Jacobs 
1980).  An  extension  of  the  TSTM  to  solid  and  fluid-filled  vertical  wall  tar¬ 
gets  has  also  been  accomplished  (Hodge  1989) . 

18.  An  important  application  of  the  temperature  prediction  models  is 
the  determination  of  the  thermal  IR  radiances  of  targets  and  backgrounds  as  a 
preliminary  step  for  the  determination  of  the  apparent  IR  scenes  that  would  be 
viewed  by  a  missile  seeker  (Bornemeier,  Bennett,  and  Horvath  1969;  Gillespie 
and  Kahle  1977;  Rounds,  Zavodny,  and  Mazzer  1980;  Zavodny  and  Mazzer  1980).  A 
particularly  useful  thermal -imaging  research  program  is  the  Air  Force  Infrared 
Modeling  and  Analysis  (IRMA)  study  that  developed  the  computer  programs 
ENVIRON/SCNGEN,  which  construct  a  target- in-background  IR  scene  as  viewed  in 
the  focal  plane  of  the  missile  seeker  optics  (Bo "kin  et  al.  1981,  Baird  et  al. 

1982) .  ENVIRON  calculates  target  and  background  temperatures,  whereas  SCNGEN 
produces  the  apparent  images  of  the  target  and  background.  The  IRMA  study  has 
been  applied  to  modeling  of  military  vehicles  such  as  tanks  (Botkin  et  al. 
1981,  Baird  et  al.  1982). 

19.  The  environmental  information  required  as  input  data  for  models 
that  describe  laser  light  and  radar  wave  scattering  from  terrain  backgrounds 
consists  of  weather  data,  terrain  type  and  land  use,  terrain  roughness,  soil 
and  vegetation  moisture  content,  terrain  roughness,  and  the  electromagnetic 
properties  of  the  terrain  background  such  as  average  dielectric  constant, 
standard  deviation  of  dielectric  constant,  and  the  average  value  and  standard 
deviation  of  the  magnetic  permeability.  The  electromagnetic  parameters  of 
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terrain  background  materials  generally  depend  strongly  on  the  terrain  moisture 
content,  which  in  turn  is  related  to  weather  and  seasonal  conditions.  The 
laser  light  bidirectional  reflectance  and  the  radar  wave  scattering  differen¬ 
tial  cross  section  can  be  calculated  from  these  parameters  (Ishimaru  1978;  Lee 
and  Kong  1985a, b). 

20.  This  report  summarizes  the  terrain  and  target  models  and  battle¬ 
field  aids  that  have  been  developed  at  WES  to  support  the  US  Army  high-value 
target  camouflage  and  target  location  and  surveillance  programs.  These  models 
have  been  created  to  satisfy  one  of  the  US  Army's  key  operational  capabili¬ 
ties:  the  detection  of  enemy  targets  located  in  complex  terrain  backgrounds. 
Reconnaissance,  Surveillance  and  Target  Acquisition  and  Automatic  Target 
Recognition  are  rapidly  being  developed  to  improve  US  Army  target  location 
capability. 


Objectives 


21.  The  objectives  of  this  study  were  to 

a.  Summarize  WES  target  and  terrain  background  signature  models. 

b.  Summarize  WES  battlefield  aids. 


Scope 


22.  The  scope  of  this  report  includes  the  following: 

a.  A  summary  of  the  thermal  signature  models  developed  to  predict 
the  physical  temperature  and  radiance  of  terrain  backgrounds. 

b.  A  summary  of  the  rough -surface  temperature  and  radiance  predic 
tion  model. 

£.  A  description  of  the  target  models  developed,  including  fluid- 
filled  targets. 

d.  The  development  of  radar  and  laser  scattering  models  for  ter¬ 
rain  backgrounds. 

e.  A  description  of  the  color -coded  regional  thermal,  laser,  and 
radar  terrain  signature  maps. 

£.  A  summary  of  the  Camouflage  Effectiveness  Evaluation  System. 

g.  A  description  of  a  thermal  target  and  terrain  scene  generator. 
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PART  II:  THERMAL  MODELS 


Two-Dimensional  Terrain  Surface  Temperature  Model 

23.  The  terrain  surface  temperature  model  TSTM  (Balick  et  al.  1981), 
developed  by  WES,  is  based  upon  a  finite-difference  solution  to  the  thermal 
diffusion  differential  equation  in  one  dimension.  TSTM  is  1-D  in  the  sense 
that  all  heat  flow  is  vertical  (e.g.,  down  into  or  up  from  the  soil)  with  no 
horizontal  effects  considered.  This  assumption  is  sufficient  in  many  cases 
where  horizontally  homogeneous  materials  are  present,  such  as  a  farmer's 
field.  Even  near  edges,  1-D  thermal  conduction  accurately  predicts  thermal 
behavior  as  long  as  one  considers  locations  at  a  distance  from  the  edge. 
However,  many  situations  require  modeling  the  horizontal  as  well  as  vertical 
heat  flow,  such  as  on  or  very  near  edges  or  for  cases  in  which  a  subsurface 
object  is  too  small  to  be  represented  accurately  by  a  1-D  calculation.  To 
overcome  these  shortcomings  of  TSTM,  a  2-D  finite  difference  thermal  model, 
named  TWOD,  was  developed.  This  model  requires  a  relatively  sparse  amount  of 
input  data,  as  in  TSTM,  but  allows  simulations  of  horizontally  discontinuous 
materials . 

Model  description 

24.  The  terrain  temperature  model  described  here  represents  an  exten¬ 
sion  to  the  1-D  TSTM  model.  The  model  calculates  a  horizontal  as  well  as 
vertical  temperature  profile  by  determining  energy  transfer  in,  out,  and 
through  the  system.  Because  the  model  is  2-D,  a  basic  assumption  is  that  the 
materials  are  uniform  in  the  third  spatial  dimension.  TWOD  predicts  a  temper¬ 
ature  profile  for  systems  that  may  be  multilayered  in  both  the  horizontal  and 
vertical  dimensions,  as  seen  in  Figure  1  with  four  materials- -A,  B,  C,  and  D. 

25.  Two-dimensional  conduction  equation.  The  temperature  profiles 
result  from  solving  the  2-D  thermal  conduction  equation  given  as 


dTix.y,  t) 
~^5t - 


a(x,y) 


d*T(x,y,  t)  +  &T(x,y,  t) 
dx2  dy2 


(1) 


where  the  observable  surface  is  y  -  0  ,  the  lower  surface  is  y  -  B  ,  the 
left  boundary  is  x  -  0  ,  and  the  right  boundary  is  x  -  R  .  The  diffusivity 
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Surface,  y  =  0 


Right 
x  =  R 


Bottom,  y  =  B 

Figure  1.  Four-material  system 

is  a(x,y)  ,  and  bit  ,  lJt  ,  rJt  ,  and  Bit  denote  heat  fluxes  at  the  appro¬ 
priate  time  t  at  the  respective  boundary  (Solomon  1982) . 

26.  Top  surface  boundary  conditions  are  calculated  in  the  same  manner 
as  in  TSTM  (Balick  et  al.  1981)  with  a  single  surface  boundary  condition  for 
all  surface  materials.  This  assumes  that  the  most  significant  heat  fluxes  at 
the  surfaces  are  vertical.  All  other  surface  boundary  conditions  are  calcu¬ 
lated  using  the  parameters  of  the  material  present  at  each  boundary  grid  point 
specification  of  one  of  three  conditions  at  the  boundary,  as  in  TSTM.  One 
should  note  that  a  combination  of  different  boundary  conditions  is  allowed  for 
the  vertical  boundaries  and  the  bottom  boundary.  A  grid  point  is  the  (x,y) 
location  of  a  numerical  value  for  the  material  temperature.  Since  the  numeri¬ 
cal  solution  is  discrete  and  not  continuous,  the  points  at  which  calculations 
are  made  define  a  2-D  grid  of  points  throughout  the  system. 

27.  System  definition  for  numerical  solution.  As  stated  earlier,  the 
model  views  the  terrain  as  a  plate  or  finite  2-D  slice  of  a  three-dimensional 
(3-D)  terrain.  In  the  case  of  a  multiple  material  in  the  vertical  and  hori¬ 
zontal  directions,  one  must  specify  the  material  properties  in  a  manner  that 
may  place  artificial  divisions  through  a  single  material.  This  is  necessary 
for  implementation  of  the  2-D  numerical  solution.  For  example,  Figure  1 
illustrates  a  system  with  four  materials,  but  for  use  with  the  model,  the 


Left 
x  =  0 


A 

B 

C 

D 
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system  must  be  specified  with  nine  blocks,  as  in  Figure  2.  One  of  the  blocks 
represents  material  A,  five  are  specified  for  material  B,  two  for  material  C, 
and  one  for  material  D.  Thermal  parameters  are  specified  for  each  block  by 
entering  the  row  and  column  number  for  each  set  of  parameters .  The  intent  of 
the  relabeling  is  to  view  the  entire  system  as  a  collection  of  rectangular 
blocks.  Care  must  be  taken  to  avoid  "thin"  blocks  that  may  violate  the  con¬ 
straints  which  ensure  numerical  stability.  This  states  that  the  product  of 
the  diffusivity  and  the  time  step  divided  by  the  sum  of  the  squares  of  the 
horizontal  and  vertical  grid  spacings  must  be  less  than  1/4  . 

Column  Number 


1 

2 

3 

1 

Row 

A 

B 

B 

Number  2 

B 

B 

B 

3 

C 

C 

D 

Figure  2.  Model  specification  for  four-material  system 


Numerical  solution  technique 

28.  Solution  within  a  block.  Within  each  block,  an  explicit  finite 
difference  scheme  is  employed  to  solve  the  2-D  heat  conduction  equation 
(Solomon  1982).  Specifically,  given  the  temperature  at  time  t,  the  new  tem¬ 
perature  at  time  t  +  At  at  the  grid  point  (x,y)  is  given  by 
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Tlx, y,  T  +  A)  *  Tlx, y,  t) 


+  a(x,y) 


(At) 


Tlx.y  +  Ay,  t)  -  2T(x,y,  t )  +  r(x,y  -  Ay,  t) 
(Ay)2 


+  T(x  +  Ax,y,  t)  -  2T(x,y,  t)  »  Tlx  -  Ax,y,  t) 

(Ax) 2 


(2) 


29.  Solution  at  an  interface.  Assuming  perfect  thermal  contact  between 
adjacent  blocks,  the  numerical  solution  equations  for  heat  transfer  across 
vertical  and  horizontal  boundaries  of  (possibly)  dissimilar  materials  can  be 
derived  by  the  principle  of  continuity  of  flux  and  temperature.  This  simply 
states  that  the  heat  flux  and  temperature  at  a  point  on  the  interface  must  be 
the  same  regardless  of  which  side  of  the  interface  the  flux  and  temperature 
are  calculated  from.  Using  this  principle,  equations  for  temperatures  at 
vertical  (Equation  3)  and  horizontal  (Equation  4)  interfaces  can  be  derived 
(Solomon  1982),  as  given  below. 


kglAxg)  ^  JCg(Ajfg) 

2aA(At)  2aB(At) 


Tlx.y,  t  *  At) 


*a(**«>  +  kB(AxB)  1 
2aA(A  t)  2aB(At)j 


Tlx.y,  t) 


Tlx  -  A x.y,  t) 


^ B 

AxB 


Tlx  +  Ax,y,  t) 


+  kg 

Axa  Axs 


Tlx,  y,  t) 


(3) 


kA^yA)  +  kg  (Ay ,) 

2aA(At)  2aB(At) 


T(x,y,  t  +  AT) 


kA(Aya)  +  JcB(AyB) 
2aA  ( A  t)  2aB(At) 


Tlx.y,  t) 


+ 


kA 
Ay a 


Tlx.y  -  Ay,  t) 


+ 


Ay* 


Tlx.y  +  Ay,  t) 


kA  +  _kg_ 
Ay *  Ay. 


T(x,y,  t ) 


(4) 
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30.  Solution  at  a  corner.  Consider  the  intersection  of  a  horizontal 
boundary  and  a  vertical  boundary,  as  in  Figure  3.  Let  T/b  denote  the  new  tem¬ 
perature  at  the  (x,y)  print  at  time  t  +  it  when  viewing  the  configuration  as  a 
vertical  interface  between  blocks  A  and  B.  Similarly,  let  Tcd  be  the  new  tem¬ 
perature  at  (x,y)  from  the  point  of  view  of  blocks  C  and  D.  Define  and  TB'D 
in  the  same  way  (Solomon  1982).  Then,  the  new  temperature  at  the  comer  point 
(x,y)  will  be 


T(x,y,  t  +  At) 


Tg*  T'cp*  t'ac  + 
4 


(5) 


Vertical 

Boundary 


u 


WES  Combined  Thermal  Model  (WESTHERM) 


31.  The  terrain  surface  thermal  models  developed  by  WES--TSTM  for  bare 
surfaces  (fialick  et  al.  1981),  VEGIE  (Balick,  Scoggins,  and  Link  1981)  for 
short  and  medium  vegetation,  and  VCTM  (Smith  et  al.  1980)  for  closed  forest 
canopies- -have  been  used  to  model  surface  and  vegetation  temperatures  for 
generation  of  surface  temperature  maps  (Weiss  and  Scoggins  1989),  among  other 
applications.  TSTM  and  VEGIE  are  traditionally  described  as  separate  models, 
even  though  they  constitute  different  sections  of  the  same  program,  because 
they  represent  two  distinct  solutions  to  the  surface  energy  budget  of  two 
types  of  system.  Each  calculates  environmental  factors  peculiar  to  the  sur¬ 
faces  for  which  it  was  designed  and  passes  the  correct  surface  temperature 
boundary  condition  to  a  single  solution  of  the  heat  conduction  equation  for 
propagation  of  the  energy  into  the  ground. 

32.  Unlike  the  TSTM-VEGIE  program,  VCTM  is  a  separate  set  of  programs 
used  to  represent  a  forest  canopy  as  a  three-layer  system,  with  additional 
layers  for  the  ground  and  sky.  One  of  the  shortcomings  of  VCTM,  at  least 
concerning  its  application  to  temperature  map  modeling  at  WES,  is  the  neces¬ 
sity  to  provide  ground  temperature  under  the  canopy  as  an  input.  This  situa¬ 
tion  occurs  because  VCTM  does  not  perform  heat  conduction  calculations  for  the 
soil  and  thus  cannot  model  the  ground  temperature  under  the  canopy.  For  ini¬ 
tial  surface  temperature  map  generation,  this  deficiency  has  been  circumvented 
by  using  VEGIE  to  model  the  ground  with  a  high  grass  cover  and  then  feeding 
the  ground  temperature  output  from  VEGIE  into  VCTM.  Such  a  temporary  solution 
provided  a  means  of  estimating  forest  canopy  temperatures  but  did  not  address 
the  true  physics  of  the  system. 

33.  To  overcome  this  limitation,  VCTM  must  be  combined  with  the  TSTM- 
VEGIE  program  in  such  a  way  as  to  allow  the  calculation  of  the  surface  energy 
budget  of  the  ground  under  the  canopy.  This  solution  requires  a  coupling  of 
the  ground  surface  temperature  to  the  canopy  temperature  in  addition  to  the 
canopy- to -ground  coupling  that  is  already  present  in  VCTM.  Following  is  a 
description  of  how  VCTM  was  combined  with  the  TSTM-VEGIE  program  to  accomplish 
this  coupling  and  thus  produce  a  combined  surface  thermal  model  (WESTHERM)  in 
which  the  canopy  is  coupled  with  the  ground. 
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Technical  approach 

34.  Background .  VCTM  models  canopy  layer  temperatures  by  performing 

energy  budget  calculations  on  each  layer  (Smith  et  al.  1980).  The  energy 

budget  calculations  consist  of  balancing  the  energy  input  to  a  layer  against 

energy  lost  by  the  layer.  Mechanisms  by  which  energy  transfers  occur  are 

Solar  radiation  (insolation)  -  Inputs  energy 
Atmospheric  longwave  radiation  -  Inputs  energy 
Canopy  layer  graybody  radiation  -  Removes  energy 
Sensible  (convection)  transfer  -  Inputs  and  removes  energy 
Latent  (evaporation)  transfer  -  Inputs  and  removes  energy 

These  same  terms  must  be  calculated  at  ground  level  under  the  canopy  to 
perform  the  ground-energy  budget.  Of  the  five,  sensible  and  latent  transfers 
are  by  far  the  most  complicated  because  of  the  nature  of  airflow  in  the  can¬ 
opy.  Fortunately,  they  are  probably  of  lower  significance  compared  with  the 
radiation  terms  to  the  undercanopy  ground  temperature.  Therefore,  these  terms 
are  calculated  using  the  VEGIE  equations  for  a  ground  surface  under  medium 
vegetation,  with  an  appropriate  reduction  of  wind  speed  under  the  canopy.  The 
VCTM  modification  work  described  here  concentrates  on  the  radiation  terms  of 
the  energy  budget  that  account  for  the  greatest  amount  of  energy  transfer. 

35.  One  of  VCTM's  major  strengths  is  its  detailed  calculation  of  radia¬ 
tion  incident  on  the  canopy  from  the  sun  and  atmosphere,  and  the  exchange 
and/or  loss  of  radiation  from  the  layers  because  of  graybody  radiation.  These 
calculations  use  three  sets  of  parameters  calculated  by  a  preprocessor  routine 
that  is  part  of  the  VCTM  package.  One  of  these  sets  is  a  matrix  W  of  proba¬ 
bilities  -of-  gap,  which  is  a  function  of  the  angle  from  vertical  for  the  forest 
canopy  used  in  the  solar- loading  calculation.  The  second  is  a  matrix  S  that 
relates  the  amount  of  thermal  energy  emitted  from  all  layers  absorbed  by  a 
layer  of  the  canopy.  The  third  set  consists  of  effective  shortwave  absorptiv- 
ities  for  each  of  the  canopy  layers.  The  absorptivities  are  calculated  using 
a  Monte  Carlo  algorithm  that  traces  the  path  of  light  from  the  source  until 
extinction.  Input  required  for  all  these  calculations  includes  (for  each 
canopy  layer)  leaf  angle  probabilities,  leaf  area  index  (LAI),  and  a  0  to  1 
factor  describing  the  degree  of  leaf  clumping.  These  values  can  be  obtained 
by  measurements  or  from  theoretical  distributions.  The  original  preprocessor 
calculated  these  values  considering  the  ground  as  a  graybody  radiation  source 
only.  However,  for  ground-energy  budget  calculations,  the  longwave  transfer 
matrix  is  modified  so  that  the  ground  is  also  a  sink  for  thermal  radiation. 
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36.  Longwave  radiation  calculation.  The  longwave  transfer  matrix  S 
calculates  the  amount  of  energy  emitted  by  each  layer  (canopy,  sky,  and 
ground)  absorbed  by  all  three  canopy  layers,  including  self-absorption  of  the 
canopy  layers.  The  elements  of  S  are  based  on  the  probability  of  gaps 
occurring  in  each  canopy  layer  and  the  mean  orientation  of  leaves  of  a  layer 
(Smith  et  al.  1980).  The  gap  probabilities  are  calculated  using  Equation  6, 
with  N(i)  being  the  LAI  of  layer  i  and  g(i,0r)  the  mean  leaf  area  layer  i 
projected  in  the  dz  direction.  This  equation  assumes  that  the  gap  probabil¬ 
ities  are  independent  of  azimuth  angle  4>z  . 


P(r, 6r)  »  <=<•> 


«/2 


gU,ez)  =  j  k{6r.6k)  fikdQk 

0 

fi*  =  -|  cos  (0k)  cos  (0r)  ,  flj,  s  |  -  ef 
Jc(0,.0*)  =  cos  (8*)  cos  ( Qz )  ^  -  tan  (♦*)  ,  0*  fc  -  8r 

=  cos"1  [-cot  (0*)  cot  (0r)] 


(6) 


37.  Changes  in  the  S  elements  are  necessary  so  that  the  ground  can  be 
a  sink  or  absorbing  layer  as  well  as  a  source  layer.  Thus,  S  will  have  five 
source  layers  (consisting  of  the  sky,  three  canopy  layers,  and  the  ground)  and 
four  sink  layers  (consisting  of  three  canopy  layers  and  the  ground) ,  instead 
of  only  the  original  three  canopy  layers.  The  same  equation  is  used  to  deter¬ 
mine  the  ground-sink  S  probability  factors,  except  that  now  the  orientation 
of  the  ground  layer  replaces  the  mean  leaf  orientation  and  the  area  projected 
is  unity.  Currently,  the  calculation  is  made  for  a  horizontal  ground,  but 
generalization  to  arbitrary  ground  slope  is  straightforward.  Including  the 
ground  in  the  matrix  is  accomplished  by  examining  the  gap  probabilities 
P(i,0r)  defined  in  Equation  7  below.  The  probability  that  light  will  pass 
through  canopy  layers  and  reach  the  ground  is  the  product  of  the  probabili- 
ties-of-gaps  for  each  layer.  Since  thermal  radiation  can  be  emitted  and 
absorbed  by  the  same  layer,  an  additional  probability  P'  is  required  which 
is  the  probability  of  a  gap  through  one-half  the  thickness  for  a  given  layer. 
This  is  given  as  (Smith  et  al.  1980) 
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Pi(r,6r)  *  Jpi  (r,Br) 


(7) 


38.  To  illustrate  the  overall  calculation,  consider  the  probability 
that  thermal  radiation  from  layer  2  will  hit  the  ground  from  a  direction  9r  . 

Pg( Qz)  =  P( 2,6r)  x  P(3,6r)  -  P'  (2,0r)  x  P( 3,6,)  x  P(2,6r)  (8) 

39.  The  total  radiation  reaching  the  ground  from  layer  2  is  then  found 
by  integrating  the  radiance  contribution  from  layer  2  over  the  hemisphere 


Lz  g  =  2x  x  fPg(Q)  x  Lj  x  |  £  •  f  |  x  d0 
o 


=  2n  x  1*2  x  V '  Pg(0)  x  Cyjj 
f^o 


(9) 


where  L2  is  the  radiation  of  canopy  layer  2  as  a  function  of  layer  tempera¬ 
ture  and  &  and  r  are  the  ground  normal  and  source  direction  unit  vectors, 
respectively.  Note  that  the  azimuthal  integral  is  factored  out  since  the 
probabilities  are  not  a  function  of  azimuth,  resulting  in  the  2ir  .  Since  L2 
is  considered  a  constant  for  a  given  layer,  it  can  also  be  brought  out  of  the 
integral,  with  the  result  being  the  S  element.  In  the  VCTM  preprocessor, 
the  integral  is  replaced  by  a  sum  over  discrete  sectors  approximating  the 
hemisphere,  and  the  parameter  CiJk  accounts  for  the  area  projection  of 
hemispheric  sector  ij  onto  layer  k.  The  total  energy  reaching  the  ground  is 
then  the  sum  of  contributions  from  this  and  the  other  two  layers  calculated  in 
the  same  manner.  This  approach  does  not  account  for  energy  emitted  from  the 
ground  that  is  reflected  back  to  the  ground,  but  since  leaves  have  a  thermal 
absorptivity  close  to  one,  multiple  reflections  should  not  introduce  a  large 
error. 

40.  Solar  loading  calculation.  In  addition  to  the  longwave  transfer 
calculation,  the  ground  energy  budget  requires  solar  loading,  or  insolation. 
Fortunately,  the  mechanism  already  exists  in  VCTM  for  this  calculation.  The 
view  angle  matrix  W  was  originally  used  to  determine  the  amount  of  radiation 
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emitted  from  the  ground  and  canopy  layers  that  exited  the  canopy  at  the  top  at 
an  angle  8  with  vertical.  This  concept  can  be  reversed  to  calculate  the 
amount  of  solar  radiation  that  reaches  the  ground.  The  calculation  requires 
only  those  elements  of  W  that  are  used  for  the  ground  contribution  to  the 
total  canopy -emitted  radiation.  Thus,  solar  contribution  to  the  ground  is 
broken  into  two  parts:  a  direct  component  and  a  diffuse  component.  Where 
available,  these  components  may  be  entered  as  input.  If  only  total  solar 
radiation  is  available,  a  number  representing  the  fraction  of  diffuse  to 
direct  must  be  estimated  (the  default  is  0.25). 

41.  Given  direct  and  diffuse  values,  total  insolation  on  the  ground 
under  the  canopy  is  given  by 


I dir  =  ^diz  x  ^(4,0<un) 


(10) 


where  Idlz  and  Sdic  are  the  insolation  and  direct  solar  components,  respec¬ 
tively,  and  W(4 , 0aun)  is  the  W  matrix  element  for  the  ground-view  factor  at 
the  solar  elevation  angle  8axm  . 

42.  The  diffuse  contribution  requires  a  summation  over  the  hemisphere 
and  is  given  by 


Idtt  m  2 


S^if  x  W(i,  6) 


(ID 


where  Idl£  and  Sdl£  are  the  diffuse  terms  and  W(4 ,8)  is  the  terra  of  W 
for  discrete  values  of  8  .  Note  that  the  summation  is  over  0  to  n/2  and  is 
multiplied  by  2  since  the  diffuse  radiation  is  considered  symmetric  about  the 
vertical  axis.  The  total  insolation  I  is  then  the  sum  Idir  and  Idi£  . 

43.  Energy  budget  for  the  ground.  Having  calculated  the  longwave  radi¬ 
ation  on  the  ground  from  the  sky  and  canopy,  L  ,  and  the  total  insolation,  I, 
the  ground  energy  budget  can  be  calculated  by  the  approach  of  TSTM  for  bare  . 
surfaces  (Balick  et  al.  1980).  This  surface  boundary  condition  is  then  propa¬ 
gated  through  the  subsurface  of  the  ground  by  invoking  the  heat  conduction 
equations  of  TSTM.  By  repeated  Invocation  of  the  modified  VCTM  module  of 
WESTHERM  to  calculate  undercanopy  ground- energy  budget  and  canopy  layer 
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temperatures ,  the  total  system  can  be  modeled  in  the  same  manner  as  the  VEGIE 
module  of  TSTM  simulates  the  effect  of  short  to  medium  vegetation  on  surface 
temperatures . 
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PART  III:  ROUGH  TERRAIN  SURFACE  TEMPERATURE  MODEL 


44.  This  section  describes  development  of  a  simple  2-D  model  for  esti¬ 
mating  the  effects  of  a  rough  terrain  surface  on  the  values  of  the  background 
radiance  that  would  be  measured  by  a  sensor  oriented  at  a  specified  depression 
angle.  This  is  done  first,  for  the  case  of  solar  loading,  by  developing  a 
two- temperature  model  (for  sunlit  and  shaded  areas)  that  predicts  the  effec¬ 
tive  brightness  temperature  of  a  rough  terrain  area.  In  addition  to  atmo¬ 
spheric  attenuation,  the  observed  background  radiance  depends  on  the  degree  of 
terrain  roughness,  texture  type,  sensor  depression  angle,  solar  altitude  angle 
(for  insolation),  and  laser  source  depression  angle  for  the  case  of  an  active 
laser  system.  Secondly,  a  similar  model  will  be  developed  to  describe  the 
ground  backscatter  from  an  active  laser  remote  sensing  system.  Both  of  these 
model  developments  are  based  on  a  2-D  calculation  in  which  the  source,  sensor, 
and  terrain  element  are  in  the  same  plane  and  only  the  depression  angles  of 
the  source  and  sensor  are  used.  Future  development  will  include  an  azimuthal 
angle  dependence  in  order  to  produce  complete  3-D  models. 

45.  The  energy  flux  measured  at  the  sensor  originates  only  from  the 
terrain  areas  accessible  to  the  sensor's  field  of  view.  Some  portion  of  a 
terrain  area  may  be  hidden  from  view  because  of  the  geometry  (amplitude  and 
wavelength)  of  the  rough  terrain  surface;  these  areas  will  not  contribute  to 
the  energy  flux  measured  at  the  sensor.  The  procedure  for  calculating  the 
backscattered  radiance  of  rough  terrain,  for  both  solar  and  laser  loadings, 
involves  determining  the  percent  of  the  viewed  terrain  surface  area  that  is 
sunlit  or  laser  illuminated  as  the  case  may  be.  For  the  case  of  insolation, 
the  sunlit  and  shadowed  portions  of  a  viewed  area  contribute  different  energy 
fluxes  because  they  have  different  surface  temperatures.  For  the  case  of  an 
active  laser  system,  the  shadowed  portions  of  the  viewed  area  make  no  contri¬ 
bution  to  the  energy  flux  measured  at  a  distant  sensor.  The  bidirectional 
reflectivity  is  interpreted  to  be  proportional  to  the  fraction  of  the  viewed 
area  illuminated  by  laser  light.  In  the  following  sections,  a  formalism  is 
developed  for  determining  the  fraction  of  a  viewed  terrain  area  that  is  illu¬ 
minated  or  shadowed. 

46.  The  model  calculations  are  performed  for  sinusoidal  terrain  sur¬ 
faces  but  can  be  extended  to  surfaces  of  arbitrary  form  by  Fourier  synthesis. 
The  macroterrain  elevations  of  real  surfaces  can  be  obtained  from  Defense 
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Mapping  Agency  databases,  the  mesoterrain  roughness  by  laser  profiling,  and 
the  microterrain  roughness  from  soil  and  vegetation  samples. 

Model  for  Solar  Illumination  on  Rouph  Terrain  Surface 

47 .  The  apparent  terrain  background  temperature  that  a  sensor  measures 
is  generally  called  the  brightness  temperature.  This  temperature  differs  from 
the  physical  temperature  of  the  ground  for  three  reasons:  atmospheric  attenu¬ 
ation,  reflected  electromagnetic  energy,  and  rough  surface  effects.  The 
reduction  in  surface  brightness  temperature  due  to  the  atmospheric  attenuation 
is  modeled  by  Electro-Optical  Systems  Atmospheric  Effects  Library.  The  effect 
of  reflected  energy  on  the  brightness  temperature  of  terrain  surfaces  has  been 
treated  in  earlier  literature.  Weiss  and  Scoggins  (1986)  presented  a  model 
for  predicting  the  effects  of  surface  roughness  on  the  brightness  temperature. 

48.  The  brightness  temperature  of  a  rough  terrain  surface  depends  on 
the  following  factors:  solar  altitude  angle,  sensor  depression  angle,  average 
terrain  roughness  amplitude,  and  average  terrain  roughness  wavelength  (Bass 
and  Fuks  1979,  Carroll  1982,  Mahrer  1982,  and  Weiss  and  Scoggins  1986).  These 
factors  will  determine  the  fraction  of  the  total  terrain  surface  area  that  is 
sunlit,  the  fraction  of  the  total  terrain  surface  area  that  is  viewed  by  the 
sensor,  and  the  fraction  of  the  viewed  area  that  is  sunlit  or  shaded.  The 
fraction  of  the  viewed  area  that  is  sunlit  and  shaded,  combined  with  the  phys¬ 
ical  temperatures  of  the  sunlit  and  shaded  areas,  is  used  to  calculate  the 
brightness  temperature  of  rough  terrain. 

49.  Initially  the  terrain  surface  is  assumed  to  have  a  sinusoidal  vari¬ 
ation  of  elevation.  The  generalization  to  surfaces  of  arbitrary  shape  can  be 
accomplished  by  Fourier  synthesis.  A  first  step  toward  calculating  the 
brightness  temperature  of  rough  terrain  is  to  determine  the  fraction  *  f  a 
sinusoidal  surface  that  is  sunlit.  This  is  done  by  referring  to  Figure  4, 
which  shows  the  intersection  of  a  sun  ray,  having  a  solar  altitude  angle  0,  , 
with  a  terrain  elevation  sinusoid.  The  intersection  points  x2,  and  x3, 
shown  in  Figure  4  can  be  used  to  calculate  the  sunlit  and  shaded  fractions  of 
the  terrain  surface  as  follows: 
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Figure  4.  Calculation  of  illumination  function 


F,  =  1  Q't  <  6,  <  -| 


(12 


F. 


L  *  *3.  -  *3, 

L 


o  <  e,  <  e; 


(13 


-  1  -  F. 


(14 


where 

F,  -  fraction  of  terrain  wavelength  that  is  sunlit 

6\  -  limiting  solar  altitude  angle  beyond  which  entire  terrain 
surface  is  sunlit 

8,  -  solar  altitude 

L  -  wavelength  of  terrain  sinusoid 

F,h  -  fraction  of  terrain  wavelength  that  is  shaded 

Therefore,  to  determine  the  fraction  of  sunlit  area  in  terms  of  the  solar 

altitude  angle,  the  terrain  roughness  amplitude,  and  the  terrain  roughness 


wavelength,  the  intersection  points  x2i  and  x3,  have  to  be  expressed  in 
terms  of  these  parameters. 

50.  The  calculation  of  the  intersection  points  x2s  and  x3s  proceeds 
as  follows.  The  terrain  surface  sinusoid  is  written  as 

y  ■  A  sin  ( kx )  (15) 

where 

A  -  amplitude  of  terrain  roughness  sinusoid 
k  -  2*/L  -  wave  number  of  terrain  roughness  sinusoid 

51.  The  slope  of  the  terrain  roughness  sinusoid  at  any  point  is  given 
by 

=  kA  cob  (kx)  (16) 

Figure  4  shows  that  if  x£„  is  defined  as  the  point  of  intersection  of  the 
tangent  sun  ray  (of  slope  tan  0a)  with  the  sinusoid,  then 

tan  de  =  kA  cos  (kx£f)  (17) 


x3M  m  L  *  X], 


(18) 


so  that 


(19) 
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52.  Equation  19  determines  the  point  x3,  .  Also  from  Equation  19  one 
concludes  that  the  limiting  solar  angle  B\  beyond  which  the  entire  surface 
is  sunlit  is  given  by 


tan  Q't  =  kA 


(20) 


53.  To  calculate  the  point  x2l  ,  one  writes  the  equation  of  the 
straight  line  representing  the  sun  ray  as  follows: 


y  =  mx  *  b 


(21) 


where  m  is  the  slope  and  b  is  the  y  intercept  at  x  -  0  (see  Figure  4) . 
The  problem  is  to  evaluate  m  and  b  in  terms  of  x3s  .  This  is  done  by 
first  noting  that 


m  =  tan  0# 


A  sin  (kx3a) 

d  +  X3, 


(22) 


where  d  is  the  positive  value  of  x  intercept  at  y  -  0  .  The  value  of  d 
can  be  obtained  from  Equation  22  to  be 

d  »  cot  sin  (kx3g)  -  x3,  (23) 


54.  Then,  from  Figure  4,  it  follows  that 

b  -  d  tain  8, 


(24) 


and  therefore, 


b  ®  A  sin  (iutj,)  -  x3#  tan  0, 


(25) 
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55.  The  equation  of  the  tangent  sun  ray  is  then  given  by  Equations  21, 
22,  and  25  as 


y  =  (x  -  x3s )  tan  0#  +  A  sin  (Jcx3#) 


(26) 


The  value  y  at  the  point  x2,  is  then  given  by 


A  sin  (kxag)  *  (x2<  -  x3<)  tan  0#  +  A  sin  (*x3-)  (27) 

56.  Equation  27  relates  x2s  and  x3,  ;  therefore,  the  simultaneous 
solution  of  Equations  19  and  27  gives  the  values  of  x2s  and  x3s  in  the 
following  functional  forms: 


x2b  =  x3f  (L#  A,0ff) 


(28) 


X3B  -  X3#(L,A,0,) 


(29) 


57.  The  simultaneous  solution  of  Equations  19  and  27  is  easily  obtained 
numerically  using  a  computer.  Then,  Equations  13  and  14  determine  the  sunlit 
fraction  F,  and  shaded  fraction  F,h  in  terms  of  L  ,  A  ,  and  6,  . 

F,  -  Fg(L, A,Qm)  (30) 

r*  -  F^U.A.Q.)  (31) 

Figures  5  and  6  show  the  dependence  of  F,  on  8,  and  L  ,  respectively. 
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Sensor  View  Model 


58.  The  construction  of  a  sensor  view  model  requires  the  development  of 
two  capabilities:  (a)  the  prediction  of  the  fraction  of  the  rough  terrain 
area  that  is  viewed  by  a  sensor  for  a  specified  depression  angle  (view  alti¬ 
tude  angle)  and  (b)  the  prediction  of  the  fraction  of  the  viewed  area  that  is 
sunlit  (or  illuminated,  in  the  case  of  laser  light).  The  fraction  of  the 
viewed  area  that  is  lighted  and  shaded  determines  the  brightness  temperature 
for  solar  loading  on  rough  terrain  and  the  bidirectional  reflectance  for  the 
case  of  laser  illumination  of  rough  terrain. 

Fraction  of  terrain  area  viewed 

59.  An  analysis  similar  to  that  done  to  calculate  the  percent  of  ter¬ 
rain  area  that  is  illuminated  can  be  used  to  calculate  the  fraction  of  a  rough 
terrain  surface  viewed  by  a  sensor.  Let  0V  be  the  view  depression  angle  and 
x2v  and  x3v  represent  intersection  points  of  a  view  line  with  the  terrain 
surface  sinusoid  as  shown  in  Figure  7.  Then,  the  viewed  and  hidden  fractions 
of  a  terrain  area  are  expressed  as 


fv  =  i  a;  <  ev  <  *  -  e;  02) 

fv  =  L  *  Xa;  0  <  ev  <  0;  (33) 

L 


^JJV  *  1  ”  Fy 


(34) 


where 

Fv  -  fraction  of  rough  terrain  surface  that  is  viewed 

-  limiting  view  depression  angle  beyond  which  entire  terrain  sur¬ 
face  is  viewed 

6V  -  view  depression  angle  (view  altitude  angle) 

Fnv  -  fraction  of  rough  terrain  surface  that  is  not  viewed  (hidden) 
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60.  The  values  of  x2v  and  x3v  are  obtained  in  a  manner  similar  to 
that  of  x2s  and  x3a  by  solving  the  following  simultaneous  transcendental 
equations 


*3v 


L  + 


tan" 

k 


'V/PI7 


-  tan2  0„ 


tan  6„ 


(35) 


A  sin(icx2v)  *  (x2Y  -  x3Y)  tan  6V  +  A  sin  ( kx3r ) 


3v' 


(36) 


61.  These  calculations  determine  x2v  and  x3v  in  terms  of  A  ,  L  , 
and  0V  for  0  <  0V  <  dY  .  i • e . , 


X^y  *  {LfA/Qy) 


(37) 


X^y  =  X^y  (L/A/Sy) 


(38) 


62.  It  is  clear  from  Equation  35  that  the  limiting  view  angle  is 
given  by 


29 


tan  0y  *=  kA 


(39) 


such  that  when  l£  <  0V  <  it  -  (£  ,  the  entire  rough  terrain  sinusoid  will  be 
viewed  (i.e.,  there  are  no  hidden  areas). 

63.  For  the  case  when  the  view  angle  is  opposite  to  the  sun  or  to  the 
laser  source  (i.e.,  when  it  -  (£  <  0V  <  it) ,  the  intersection  points  of  the  view 
line  with  the  terrain  surface  sinusoid  are  designated  by  x4v  and  Xjv  .  The 
values  of  x*v  and  x5v  are  given  by 


(40) 


(41) 


64.  The  points  x*v  and  x5v  are  related  in  this  way  only  for  conju¬ 
gate  view  lines,  i.e.,  0V  and  it  -  0V  .  For  the  conjugate  view  lines,  the 
view  fraction  is  given  by 


Fv 


L  *  *4v  ~ 

L 


X  -  0y  <  0V  <  x 


(42) 


-  L  *  X*v  ~  X>r 
L 


0  <  0V  <  7i  -  0* 


(43) 


65.  Thus,  Fv  is  symmetrical  under  0V  -»  tt  -  9V  .  Combining  Equations 
33-36  and  Equation  42  gives  the  viewed  and  nonviewed  fractions  as 

Fv*F(A.L.bv)  (44) 
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Fov  =  Fdv(A,  L,  By) 


(45) 


for  0  <  0V  <  it  .  The  function  Fv  is  the  same  as  Fs  that  appears  in 
Figures  5  and  6. 

Fraction  of  viewed  area  illuminated 

66.  The  important  physical  quantity  to  calculate  for  evaluating  sensor 
performance  is  the  fraction  of  the  viewed  area  that  is  illuminated,  because  it 
is  this  fraction  that  determines  the  effective  brightness  temperature  of  rough 
terrain  that  a  sensor  will  detect  in  its  field  of  view,  and  it  is  this  frac¬ 
tion  that  determines  the  bidirectional  reflectance  for  the  case  of  a  laser 
source  and  detector. 

67.  Let  fvg  be  the  fraction  of  the  viewed  area  that  is  illuminated. 
Since  fv,  refers  both  to  the  source  and  the  viewer,  it  will  in  general  be  a 
function  of  both  0a  and  0V  ,  and  can  be  written  in  the  following  functional 
form: 


f„  =  f,(A.L,  8,,0V) 


(46) 


The  fraction  fva  will  be  calculated  in  terms  of  the  intersection  points 
x2»  »  x3«  •  x2v  .  x3v  .  and  also  in  terms  of  x*v  and  Xsv  for  the  case 
n  -  <  0V  <  n  .  This  will  give  the  desired  result  because  these  intersec¬ 

tion  points  have  already  been  expressed  in  terms  of  A  ,  L  ,  0a  ,  and  0V  . 

68.  Figure  8  indicates  that  the  calculation  of  must  be  separated 

into  three  view  angle  classes;  subsequently,  each  class  must  be  separated  into 
several  subclasses.  This  is  done  as  follows: 

Class  I:  0  <  0V  <  0^ 


Subclass  Ai  0  <  e„  <  6,  f v,  =  1 


(47) 


Subclass  B:  8.  <  0  <  0*  f _  = 

•  v  w  L  +  *av  -  *3V 


(48) 
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Figure  8.  Calculation  of  fraction  of  viewed 
area  that  is  illuminated 
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Subclass  A:  0j  <  6,  f „  *  1 


(49) 


Subclass  B:  0  <  8,  < 


f  =  L  *  **•  ~ 


(50) 


Class  III:  x  - 


C  <  <  * 


Subclass  A:  x2t  <  xAr  fv 


b  *  *2.  -  *3# 

Z.  +  X^y  ~  A^y 


(51) 


Subclass  B:  x3t  <  x4. 


f  =  L  *  x«v  - 
**  b  *  xtr-  xsr 


(52) 


Subclass  C :  x2#  <  x5r 


^  *  jVr  ~  ^sr  *  xu  *  x3< 

b  *  X4r  - 


(53) 
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69.  The  function  fva  is  not  symmetrical  under  the  transformation 
8V  -*  it  -  tfv  .  The  fraction  of  the  viewed  area  that  is  not  illuminated  is 
given  by 


1  -  f „ 


(54) 


Figure  9  shows  the  dependence  of  fv,  on  0V  and  0,  . 


Figure  9.  Angular  dependence  of  fraction  of 
viewed  area  that  is  illuminated 

Effective  brightness  temperature 

70.  The  effective  brightness  temperature  of  rough  terrain  as  observed 
in  the  field  of  view  of  an  infrared  sensor  can  be  calculated  in  terms  of  the 
fraction  of  the  viewed  area  that  is  sunlit,  fyg  .  Excluding  reflected  energy, 
the  radiance  of  the  sunlit  and  shaded  portions  of  a  rough  terrain  area  are 
given  respectively  by 


= 


e.C^AX 


Xs  [exp  <<vur.)  -  1] 


(55) 
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_ 'nt^AX 

X5  [exp  (C2/ (XT^)  -  1] 


(56) 


*«.  = 


where 


N, 

«. 

Cl 

AA 

A 

C2 

T, 


£ 


ns 


radiance  of  sunlit  area,  W  m”2  sr”1 
emissivity  of  sunlit  surface 
1.192  x  10®  W  pm*  m'2  sr”1 
spectral  wavelength  interval,  pm 
weighted  average  of  wavelength,  pm 
1.434  X  10*  pm  ®K 
temperature  of  sunlit  area,  *K 
radiance  of  shaded  area,  W  m”2  sr”1 
emissivity  of  shaded  surface 
temperature  of  shaded  area,  ®K 


71.  The  average  radiance  of  a  rough  terrain  surface  in  the  field  of 
view  of  a  sensor  is  given  by 


N=  fy,N,  *  (1  -  fy.)^ 


(57) 


72.  Then,  the  effective  brightness  temperature  is  defined  by 


N  = 


_ q  AA _ 

A5  [exp  ( Ca/XT0ft )  -  1] 


(58) 


where  T,ff  is  the  effective  brightness  temperature.  In  this  way  the  rough¬ 
ness  model  determines  the  following  functional  dependence: 

T,ff  *  Ttff  (59) 
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73.  The  terrain  amplitude  and  wavelength,  A  and  L  ,  respectively, 
can  be  estimated  from  statistical  terrain  roughness  parameters  (Weiss  1981, 
1982).  Figure  10  shows  an  example  of  the  dependence  of  Teff  on  the  view 
angle  0V  . 


A/L  -  1.0 


Figure  10.  Angular  dependence  of  brightness 
temperature 
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PART  IV:  INFRARED  TARGET  AND  SLOPED  TERRAIN  MODEL 


74.  This  pact  of  the  report  describes  the  procedure  for  calculating  the 
intrinsic  radiance  of  the  target  and  background  surface  facets.  The  procedure 
involves  the  calculation  of  surface  temperatures  by  using  nn  energy  budget 
condition  (Budyko  1974) .  The  radiative  energy  input  to  this  budget  for  a 
horizontal  surface  consists  of  (a)  a  short-wave  (SW)  solar  component  and  (b)  a 
direct  longwave  (LW)  atmospheric  component  (Balick  et  al.  1981).  For  vertical 
target  surfaces,  there  are  three  additional  components:  (a)  a  SW  component 
reflected  from  the  ground  to  the  target,  (b)  a  LW  atmospheric  component 
reflected  from  the  ground  to  the  target,  and  (c)  a  LW  irradiance  on  the  target 
resulting  from  the  thermal  radiance  of  the  ground.  Thus,  for  a  horizontal 
surface  there  are  two  primary  radiative  power  inputs,  whereas  for  a  vertical 
surface  there  are  five.  The  original  TSTM  was  developed  for  horizontal  sur¬ 
faces.  To  use  the  TSTM  for  predicting  the  temperature  of  a  vertical  target 
surface,  one  must  consider  the  necessary  radiative  processes  associated  with  a 
vertical  wall  (Weiss  and  Scoggins  1989) . 

75.  Because  the  IR  target  and  background  radiances  are  generally  mea¬ 
sured  in  a  specific  wavelength  interval  such  as  AX  -  3-5  or  7-14  pm  ,  the 
predicted  radiances  must  be  calculated  for  these  spectral  windows.  However, 
the  calculation  of  the  equilibrium  target  and  background  surface  temperatures 
requires  broadband  radiance  values .  Therefore ,  both  broadband  and  spectral 
window  radiances  are  calculated  in  this  report. 

Irradiance  on  Target  and  Background 

Basic  photometric  definitions 

76.  The  spectral  radiance  of  a  target  surface  is  defined  as  the  amount 
of  power  radiated  per  unit  area  per  unit  solid  angle  and  per  unit  wavelength 
as  follows  (Holter  et  al.  1962;  Kruse,  McGlaughlin,  and  McQuistan  1962; 
Jamieson  et  al.  1963): 


dP  -  Nx  cos  0  dA  dQ  dX 


(60) 


36 


where 


P  -  emitted  power,  W 

NA  -  spectral  radiance,  W  m'2  pm'1  sr'1 

6  -  angle  between  direction  of  observation  and  the  normal  to  target 
surface ,  rad 

A  -  area  of  target  surface,  m2 
Q  -  solid  angle,  sr 
A  -  wavelength,  pm 

77.  The  radiance  of  a  surface  is  defined  as 

N  =  f  Nj.dX,  W  m'2  sr'1  (61) 


and  is  associated  with  specified  wavelength  intervals.  The  target  facet 
radiances  used  in  TSTM  are  associated  with  the  3-  to  5 -pm  and  the  7-  to  14-pm 
wavebands . 

78.  The  magnitude  of  the  IR  radiance  associated  with  the  background  and 
target  is  due  to  a  number  of  physical  processes  associated  with  the  absorp¬ 
tion,  emission,  and  reflection  of  electromagnetic  energy.  The  total  IR  radi¬ 
ance  associated  with  target  and  background  surfaces  consists  of  a  thermal 
component  and  a  reflected  component.  The  thermal  radiance  depends  on  the 
absolute  temperature  of  the  target  and  background  facets. 

79.  The  irradiance  H  (W  m'2)  is  the  electromagnetic  power  that  falls 
on  a  unit  area  of  a  surface.  The  radiance,  or  emissive  power,  of  one  surface 
may  impinge  on  another,  thereby  becoming  an  irradiance  for  the  second  surface. 
It  is  important  to  distinguish  between  radiance  and  irradiance  in  the  calcula¬ 
tions  shown  in  this  part  of  the  report. 

Irradiance  on  the  ground 

80.  The  irradiance  on  the  ground  surface  consists  of  the  following  two 

terms : 

fi.  Direct  SW  insolation  on  ground  surface.  The  SW  solar  irradi¬ 
ance  on  the  ground  is  given  by  (Balick  et  al.  1981) 

-  S0Wl  cos  Z,  (62) 
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where 


i$(SW)  -  solar  irradlance  on  ground,  W  m-2 

SD  -  solar  constant  -  1,395.0  W  nf2 

Ux  -  atmospheric  transmission  factor  (Balick  et  al. 

1981) 

Z,  -  solar  zenith  angle 

1*.  Direct  LW  atmospheric  irradlance  on  ground.  The  broadband 

atmospheric  LW  irradlance  on  the  ground  is  given  by  the  Brunt 
equation  (Balick  et  al.  1981). 

Hfa,(LW)  *  ot£  (c  +  b  (63) 

where 

“  broadband  atmospheric  LW  irradlance  on  ground, 

W  m'2 

a  -  5.673  x  10'8  W  m*2  °K'* 

Ta  —  ambient  air  temperature,  °K 
c,d  -  empirical  coefficients  (defined  below) 
ea(Ta,RH)  -  water  vapor  pressure  (Balick  et  al.  1981) 

RH  -  relative  humidity 

81.  The  LW  irradlance  is  diffuse.  The  Hftja(LW)  function  contains 
empirical  coefficients  that  are  chosen  to  have  the  values  c  -  0.73  and 

b  -  0.06  .  These  values  are  slightly  larger  than  the  values  c  -  0.61  and 
b  -  0.05  given  in  Balick  et  al.  (1981).  The  new  values  of  c  and  b  are 
chosen  because  they  significantly  improve  the  output  of  the  TSTM  when  it  is 
applied  to  various  soil  types. 

82.  The  Brunt  equation  is  broadband  and  thus  covers  all  frequencies. 
The  corresponding  irradlance  on  the  ground  surface  in  the  wavelength  interval 
AA  is  given  by 


A  HZa.iLVf) 


%C f  AA 


(c  +  bJe2 

v  ^  "  (64) 


(64) 
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where 


Al^aCLW)  -  spectral  window  LW  irradiance  on  ground,  W  m'2 
Cf  -  1.192  x  10*  W  pm*  m'2  sr’x 

AA  -  wavelength  interval  -  7  and  2  pm,  respectively 
Cf  -  1.4338  x  10*  pm  °K 

Aav  -  weighted  average  wavelength;  in  the  spectral  window  7  to 
14  pm  -  10.76  pm,  and  for  3  to  5  pm  -  4.8  pm 

Irradiance  on  a  vertical  target  facet 

83.  The  TSTM  was  originally  designed  to  predict  the  ground  surface 
temperature  or  the  temperature  of  a  horizontal  surface.  This  model  has  been 
modified  to  describe  vertical  surfaces  by  including  the  important  radiative 
interchanges  between  the  vertical  target,  the  atmosphere,  and  the  background. 

84.  The  radiative  transfer  at  the  target  surfaces  has  LW  and  SW  compo¬ 
nents  and  includes  five  major  irradiance  components,  as  described  in  the  fol¬ 
lowing  paragraphs. 

85.  Direct  SW  insolation  on  target.  The  direct  SW  solar  irradiance  on 
a  target  surface  is  given  by 


tfQtar(SW)  =  S0WX  cos 


(65) 


where 

1^“(SW)  -  solar  irradiance  on  target  surface,  W  m-2 

<f>t  -  angle  between  surface  normal  and  direction  of  the  sun 

86.  A  simple  analysis  shows  that  (Balick  et  al.  1981) 


cos  <f>t  -  sin  Z8  sin  (SL)  cos  (S^  -  SL^)  +  cos  Z8  cos  (SL)  (66) 


where 

SL-  slope  angle  of  target  surface 
Saz  “  solar  azimuth  angle 
SLaz  -  azimuth  angle  of  normal  to  target  surface 

87.  The  result  in  Equation  66  follows  from  the  fact  that  the  unit 
vector  normal  to  the  surface  of  the  target  is 
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Rt  *  1  sin  ( SL )  cos  (SLju)  *  3  sin  (SL)  sin  (SLM)  +  k  cos  (SL) 


(67) 


and  the  unit  vector  in  the  direction  of  the  sun  is  given  by 

Ra  =  l  sin  Za  cos  SM  *  3  sin  Za  sin  +  k  cos  Z,  (68) 

88.  By  the  definition  of  the  scalar  product,  it  follows  that 

cos  ■  Rt  •  Rg 

«  sin  (SL)  sin  Zt  [cos  (SLj^)  cos  &AZ 

(69) 

+  sin  (SLM)  sin  S^]  ♦  cos  (SL)  cos  Zm 
=  sin  Zg  sin  (SL)  cos  (S^  -  SLt +  cos  ZB  cos  (SL) 


This  is  the  same  result  as  in  Equation  66. 

89.  SW  solar  radiation  reflected  from  background  to  target.  The  SW 
irradiance  on  the  target  surface  due  to  specularly  reflected  insolation  from 
the  ground  is  given  by 


»rfc.“  ,(S*0  *  r6(SW)SQW1  cos2  Zg  cos 


(70) 


where 

rg(SW)  -  SW  reflected  coefficient  of  ground  -  1  -  cg(SW) 
eg(SW)  -  SW  absorptivity  of  ground 

4%  -  angle  between  surface  normal  and  the  reflected  sun  ray 

90.  For  diffuse  reflection  of  solar  insolation  from  the  ground,  the  SW 
irradiance  on  the  target  surface  is  given  by 


40 


=  A  rg(SW)  S&G  COS  z, 


(71) 


where  G  Is  the  radiation  configuration  factor  that  will  be  defined  later. 

If  the  ground  is  taken  to  be  infinite  half -plane,  then  G  -  ir/2  . 

91.  Simple  geometry  shows  that 

cos  *  sin  zt  sin  (SL)  cos  (S«  ~  SLju )  -  cos  Za  cos  (SL)  (72) 


The  result  in  Equation  72  is  obtained  by  calculating  the  scalar  product 
between  the  unit  vector  normal  to  the  target  surface  and  the  unit  vector  in 
the  (negative)  direction  of  the  reflected  sun  ray.  The  unit  vector  in  the 
(negative)  direction  of  the  reflected  sun  ray  is  given  by 

Ra  *  i  sin  Zm  cos  S^,  +  j  sin  Za  sin  ~  *  cos  Za  (73) 


92.  By  definition  it  follows  that 

cos  =  iit  •  tia 

*  sin  ( SL )  sin  Za  [cos  (SL^)  cos 

(74) 

+  sin  (SLM)  sin  S^]  -  cos  Za  cos  ( SL ) 

=  sin  (SL)  sin  Za  cos  (SM  -  SL IJiZ)  -  cos  Za  cos  (SL) 


This  is  just  the  result  given  in  Equation  72. 

93.  Direct  LW  atmospheric  irradiance  on  target  surface.  The  broadband 
atmospheric  IR  irradiance  on  a  horizontal  target  surface  is  assumed  to  be 
given  by  the  Brunt  equation  as  in  the  case  of  the  LW  irradiance  on  the  ground 
(Equation  63). 
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h£UW)  »  Hftm ( AW) 


(75) 


where  AH,ta(LV)  is  the  direct  broadband  LW  irradiance  on  target  surface,  in 
W  in"2. 

94.  For  a  vertical  target  surface,  the  broadband  atmospheric  IR  irradi 
ance  is  given  by 


h££(LW)  =  ~  J&.OJV) 


(76) 


The  corresponding  atmospheric  irradiance  on  the  horizontal  target  surface  in  a 
spectral  range  is  assumed  to  be  the  same  as  spectral  atmospheric  irradiance  on 
the  ground,  as  given  by  Equation  64,  so  that 

AH^(LW)  =  AH*m(LW)  (77) 

where  AH^(LW)  is  the  spectral  window  LW  atmospheric  irradiance  on  the  tar¬ 
get,  in  W  m-2. 

95.  For  a  vertical  target  surface,  the  atmospheric  spectral  irradiance 
is  given  by 


AH^(LW)  =  1  AH^ta(LW) 


(78) 


96 .  LW  atmospheric  radiation  reflected  from  ground  to  target.  Because 
the  atmospheric  LW  radiation  is  diffuse  in  nature,  the  radiation  reflected 
from  the  ground  will  also  be  diffuse,  and  the  LW  broadband  irradiance  on  a 
raised  target  surface  is  given  by  (General  Electric  Company  1975) 


i  rrNf'tiw  cos  et  cos  eff  JB  J. 

JJ - p - 2 


(79) 


where  the  diffuse  reflected  LW  ground  radiance  is  given  by 
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HftmlLM 


(80) 


N*0i(LW) 


rg(LW) 

K 


and  where  h&n(LW)  is  given  by  Equation  63  and 
At  -  target  area 

6t  -  angle  between  normal  to  target  and  reflected  ray 

6 g  -  angle  between  normal  to  ground  and  incident  ray 

r  -  variable  distance  between  points  on  ground  and  on  target 
surface 

Ag  -  area  of  ground 

rg(LW)  -  LW  reflectively  of  ground  -  1  -  eg(LW) 

«g(LW)  -  LW  absorptivity  of  ground  -  LW  emissivity  of  ground 

97.  For  the  case  when  the  thermal  ground  radiance  ^er(LW)  is  inde¬ 
pendent  of  the  position  in  the  ground  plane,  the  integral  in  Equation  79  can 
be  evaluated  analytically  with  the  result  that,  for  a  half-plane, 

H?*f'  g  ( LW)  =  GN?.e(LW)  (81) 

where  (General  Electric  Company  1975) 


G  -  tern-1  (£>/c) 


(82) 


and 


a,b  -  length  and  width,  respectively,  of  quarter-plane 

c  -  vertical  height  of  observation  point  on  vertical  target  wall 

98.  For  the  case  of  a  homogeneous  semi- infinite  half-plane  background, 
G  -  »r/2  ,  Equation  81  becomes 


(LW) 


(83) 


99.  The  corresponding  reflected  LW  atmospheric  radiation  in  the  7-  to 
14-/jm  waveband  is  obtained  from  Equations  80  and  81  to  be 
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AH%;ig(LW)  -  G&N?,f(LW) 


(84) 


where 


AN&iLh 0 


rgUW) 

x 


AHSLlUf) 


(85) 


and  AHft^LW)  Is  given  by  Equation  64. 

100 .  LW  irradiance  on  target  due  to  thermal  radiance  of  the  ground.  The 
diffuse  LW  irradiance  on  a  raised  target  due  to  the  thermal  radiance  of  the 
ground  is  given  by  (General  Electric  Company  1975) 


^&«r(l>W)  cos  0t  cos  6, 


£  dAg  dAt 


(86) 


where  the  broadband  thermal  ground  radiance  is  obtained  from  Planck's  law  as 


NtherlLW) 


=  egum 


a  7 g 


(87) 


and 

Cf  -  1.192  x  108  W  /an*  nf2  sr'1 
q?  -  1.4338  x  10*  nm  °K 
a  -  5.673  x  10’8  W  nf2  °K'* 

Tg  -  ground  temperature,  °K 

101.  For  the  case  of  a  homogeneous  half -plane  background,  the  integral 
in  Equation  86  simplifies  to 
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H£z,gum  -  gn^alvT) 


(88) 


where  the  geometry  factor  G  is  given  by  Equation  82. 

102.  The  corresponding  spectral  radiance  of  the  ground  is  given  by 
Equation  86  with  H5j#t(LW)  replaced  by 


A  NttoALW 


eg(LW)C?  X? 


€q(LW)  Ci  LX 


(  H  > 

Is 

exp 

ca 

-  1 

Atv 

1  T 
^  A»vAgJ 

(89) 


103.  For  the  case  of  a  homogeneous  half-plane  background,  the  spectral 
irradiance  on  the  target  due  to  the  thermal  ground  radiance  is  obtained  from 
Equation  88  to  be 


Mi^T.g(LW)  «  GAn£,»t  (LW) 


(90) 


104.  Less  important  radiative  transfer  effects  are 

g.  Shortwave  solar  radiation  reflected  from  target  surface  to  the 
background  and  reradiated  as  IR  thermal  emittance,  which 
irradiates  the  target. 

t).  Atmospheric  IR  radiation  reflected  from  target  surface  to  the 
background  and  reradiated  as  a  thermal  IR  emittance,  which  in 
turn  irradiates  the  target. 

105.  The  total  SW  irradiance  on  the  target  surface  is  obtained  from 
Equations  65  and  70  to  be 


(SW)  -  {SW)  +  Hrtlg(SW) 


(91) 


The  total  LW  irradiance  on  the  target  surface  is  obtained  from  Equations  75, 
79,  and  86  to  be 
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H***  (LW)  «  h££(LW)  *  H£t.g(LW)  *  HaXr.glLm 


(92) 


106.  For  horizontal  surfaces  such  as  the  ground,  the  target  irradiances 
given  in  Equations  91  and  92  reduce  to  those  given  in  Equations  62  and  63, 
respectively. 

107.  The  absorbed  part  of  the  total  irradiance  on  a  target  surface 
enters  the  calculation  of  the  target  surface  temperature.  The  absorbed  power 
per  unit  area  is  obtained  from  Equations  91  and  92  as 

*  €t(SW)  (SW)  +  et(LW)  Hm  (LW)  (93) 


where 

H£g5  -  power  absorbed  per  unit  area,  W  m'z 
ct(SW)  «  SW  absorptivity  of  target  surface  -  SW  emissivity 
«t(LW)  -  LW  absorptivity  of  target  surface  -  LW  emissivity 

108.  It  is  this  total  absorbed  power  that  enters  the  calculation  of 
target  facet  temperature  using  the  computer  program  TSTM.  The  procedure  for 
calculating  the  target  surface  temperature  using  Equation  93  is  essentially 
the  same  as  used  for  calculating  the  ground  temperature. 

109.  For  surfaces  that  are  not  horizontal,  the  orientation  of  the  sur¬ 
face  with  respect  to  the  direct  and  reflected  radiation  components  enters  the 
irradiance  calculation.  This  introduces  a  dependence  of  the  predicted  target 
surface  temperature  on  the  slope  and  azimuth  angles  of  the  surface  and  a 
dependence  on  the  solar  zenith  and  solar  azimuth  angles. 

Prediction  of  Target  and  Background  Temperatures 

110.  The  equilibrium  surface  temperature  is  determined  by  satisfying  an 
energy  flow  budget  equation  at  the  surface  of  a  target  or  background  (Budyko 
1974,  Balick  et  al.  1981).  The  power  input  includes  the  absorption  of  a  SW 
insolation  component  and  a  LW  atmospheric  term  as  given  by  Equation  93. 

Energy  flow  budget 

111.  The  surface  temperature  of  the  target  and  background  is  calculated 
using  the  power  budget  of  the  TSTM  computer  program  (Balick  et  al.  1981). 
Several  physical  processes  are  included  in  this  model: 
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£.  Radiative  absorption  at  the  surface  (Equation  93) . 

b .  Thermal  emission  from  the  surface. 

c.  Conduction  of  heat  in  the  solid  background  or  solid  target. 

d.  Conduction  and  convection  by  the  atmosphere. 

£.  Latent  heat  of  evaporation. 

The  energy  flow  budget  for  solid  targets  and  backgrounds  has  been  fully 
described  (Balick  et  al.  1981). 

Convection  and  conduc- 
tion  in  fluids  and  solids 

112.  In  the  energy  budget  of  the  original  TSTM,  heat  conduction  is  the 
mode  of  energy  transfer  into  a  solid  target  or  background.  However,  heat 
convection  is  the  dominant  heat  transfer  process  in  fluid- filled  targets  such 
as  petroleum  storage  tanks.  Therefore,  in  order  to  use  the  TSTM  computer 
program  for  petroleum  storage  tanks,  the  actual  heat  convection  process  must 
be  converted  to  an  equivalent  heat  conduction  problem. 

113.  The  heat  conduction  process  is  described  by  two  parameters:  heat 
conductivity  and  thermal  diffusivity.  Thermal  diffusivity  is  defined  as  fol¬ 
lows  (Hodge  1989) : 


a 


(94) 


where 

a  -  thermal  diffusivity,  cm2  min"1 
k  -  thermal  conductivity,  cal  (cm  min  #K)"1 
p  -  mass  density  of  solid,  g  mass  cm'3 
Cp  -  specific  heat  at  constant  pressure,  cal  (g  °K)_1 

114.  The  heat  conduction  equation  in  one  dimension  is  (Hodge  1989) 


dT  _ 
‘ 


a 


<pt 

ax2 


(95) 


where 

T  -  temperature,  "K 
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t  -  time 

X  "  distance,  cm 

115.  Equation  95  is  a  diffusion  equation  that  describes  the  transfer  of 
heat  by  the  spreading  of  molecular  agitation  in  a  lattice.  This  description 
is  adequate  for  solids  but  not  for  fluids  like  gases  and  liquids  whose  mole¬ 
cules  are  free  to  move  in  bulk  flows . 

116.  For  fluids,  convective  heat  transfer  totally  dominates  the  heat 
conduction  process.  Equation  95  and  its  solutions  are  not  applicable  to  heat 
convective  processes.  In  fact,  the  convective  heat  transfer  equation  is 


dT 

dt 


+ 


u 


&T 

dX2 


(96) 


where  u  -  u(x,t)  -  fluid  speed  .  The  solution  of  Equation  96,  the  continuity 
equation,  and  the  Navier- Stokes  equations  yields  the  temperature  and  fluid 
velocity  distribution  in  the  fluid.  The  solution  of  these  three  equations  is 
difficult  for  realistic  fluids,  and  a  simpler  engineering  solution  to  the 
convection  problem  that  utilizes  Equation  95  with  an  effective  conductivity  is 
applied  (Hodge  1989) . 


Target  Facet  Radiance 

117.  The  total  IR  emittance  from  a  target  or  background  surface  in  a 
wavelength  interval  AA  is  the  sum  of  a  reflected  IR  component  and  a  thermal 
IR  emittance  calculated  from  the  predicted  target  surface  temperature. 

Thermal  IB- rgdlance  of  target 

118.  The  thermal  radiance  of  a  target  or  background  surface  in  a  wave¬ 
length  interval  AA  is  calculated  from  Planck’s  blackbody  radiation  law  as 
follows 
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AN^(LW)  =  €  t  { LW)  C" 


€t  (LW)  C*  AX 


(97) 


where 

ANthJr  -  thermal  radiance  of  target,  W  m-2  sr_1 
Tt  -  target  surface  temperature,  °K 

119.  The  atmospheric  IR  transmission  windows  are  3  to  5  pm  and  7  to 

14  pm.  The  important  parameters  for  the  determination  of  thermal  radiance  are 
the  target  surface  temperature,  the  emissivity  of  the  surface,  and  the  trans¬ 
mission  wavebands. 

Reflected  IR  .remand? 

120.  In  addition  to  the  thermal  IR  radiance  of  a  surface,  there  is  a 
reflected  IR  radiance  in  a  wavelength  interval  AX  ,  which  for  a  vertical  sur¬ 
face  is  due  in  part  to  the  IR  atmospheric  radiation  reflected  directly  from 
the  target  and  in  part  to  the  IR  atmospheric  radiation  reflected  from  the 
background  and  then  reflected  again  from  the  target  surface.  There  is  also  an 
IR  reflected  component  due  to  the  thermal  emittance  of  the  background,  which 
is  reflected  by  the  target  surface.  The  total  reflected  radiance  in  a  wave¬ 
length  interval  AX  is  calculated  from  the  total  irradiance  on  the  target 
surface  in  the  wavelength  interval  AX  ,  which  is  obtained  from  Equations  77, 
84,  and  90  to  be 


AW“r(LK0  -  A tf.'SJaW)  +  AH%ig{Ltn  +A 


(98) 


121.  The  reflected  IU  radiance  in  the  wavelength  interval  AX  is  cal¬ 
culated,  under  the  assumption  that  it  is  diffuse,  in  the  following  manner: 


49 


A nZUlW)  =  AH“rUW) 


where  rt(LW)  -  LW  reflectivity  -  1  -  «t(LW)  . 

Total  target  and  ground  radiance 

122.  The  total  target  facet  IR  radiance  in  the  wavelength  interval  AA 
is  obtained  finally  as  the  sum  of  the  thermal  radiance  and  the  LW  reflected 
radiance  as  follows: 


-  A NoZt(LW)  +  AA T^f(LW) 


(100) 


where  ANtgJr(LW)  is  given  by  Equation  97.  This  is  the  intrinsic  target 
radiance  that  serves  as  input  for  the  determination  of  the  apparent  radiance 
seen  by  the  IR  detector. 

123.  The  ground  radiance  can  be  obtained  from  Equations  98  through  100 
as  a  special  case 


A *  AHfu,(LW) 


(101) 


AA#/(LW)  =  rZ--L—  UH9(LV/) 


(102) 


A N*(LW)  -  A N^Z(LW)  +  A N?,f(LW) 


(103) 


where  rg(LW)  -  1  -  es(LW)  and  where  Ai^(LW)  and  Al^,.r(LW)  are  given  by 
Equations  64  and  89,  respectively. 


124.  A  blackbody  brightness  temperature  for  a  target  such  as  a  POL  tank 
can  be  associated  with  a  radiance  value  as  follows: 
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(104) 


»pB 


X-Vln 


1  + 


C*  AA ' 
AW  A*Vj 


where 

Tb  -  brightness  temperature,  °K 

C?  -  1.4338  X  10*  pm  °K 

C?  -  1.192  x  10®  W  pm4  m'2  sr*1 

AN  -  average  apparent  radiance  of  target,  W  m"2  sr"1 
AA  -  wavelength  window,  pm 

Aav  "  weighted  average  wavelength  in  window  AA 
In  this  way,  a  temperature  can  be  associated  with  each  of  the  average 
radiances  that  have  been  defined  for  a  target. 
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PART  V:  THERMAL  MODELING  OF  PETROLEUM,  OIL,  AND  LUBRICANTS 

BULK  STORAGE  TANKS 

Background 

125.  Reliable  models  are  available  that  can  predict  the  thermal  surface 
conditions  for  grass-  and  vegetation- covered  surfaces  (Balick,  Scoggins,  and 
Link  1981)  as  well  as  for  bare  soil  and  concrete  (Balick  et  al.  1981).  Com¬ 
plementary  thermal  surface-condition  models  for  fixed-installation  targets 
such  as  POL  bulk  storage  tanks  are  needed  but  are  not  generally  available. 
Because  POL  storage  tanks  are  high- asset  targets  and  because  modeling  method¬ 
ology  appropriate  for  these  tanks  could  also  be  used  for  other  fixed- 
installation  targets,  POL  tanks  were  chosen  as  the  target  to  be  thermally 
modeled. 

126.  The  objective  of  this  effort  was  the  preliminary  development  and 
validation  of  a  thermal  model  to  realistically  predict  the  surface  temperature 
histories  of  POL  bulk  storage  tanks.  To  possess  sufficient  fidelity,  such  a 
model  must  account  for  external  factors  such  as  air  temperature,  solar  insola¬ 
tion,  wind,  and  cloud  conditions  as  well  as  internal  convection  and  storage 
effects.  As  with  the  Terrain  Surface  Temperature  Model  (Balick  et  al.  1981), 
a  premium  was  placed  on  simplicity  and  flexibility.  A  model  with  these  char¬ 
acteristics  must  consider  and  emulate  the  dominant  physical  phenomena  that 
influence  the  surface  temperatures. 

127.  A  simple  but  physically  realistic  model  to  predict  the  temperature 
time -dependence  of  the  external  surfaces  on  POL  bulk  storage  tanks  is  summa¬ 
rized  in  this  section  (Hodge  1989).  Model  development  is  first  considered  in 
terms  of  the  physical  phenomena  that  influence  the  temperatures  of  each  of  the 
surfaces,  and  from  these  considerations  the  physical  bases  of  the  overall 
model  are  selected.  The  WES  vertical  wall  target  model  was  used  to  calculate 
the  radiant  energy  flux  on  the  target  wall  (Weiss  and  Scoggins  1989) . 

128.  A  lumped-mass  surface- temperature  thermal  model  for  POL  bulk 
storage  tanks  is  developed.  The  model  uses  12  ordinary,  coupled,  first-order 
differential  equations,  obtained  from  energy  balances,  to  describe  the  time- 
dependent  variation  of  the  surface  temperatures  of  a  partially  filled  POL 
tank.  Convective  heat  transfer  and  direct  and  diffuse  irradiation  are  modeled 
as  energy  transfers  into  and  out  of  the  12  lumped  masses.  One -dimensional 
transient  conduction  is  considered  in  the  soil  below  the  tank.  The  model  is 
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implemented  in  a  computer  simulation.  Comparisons  of  the  predicted  tank  sur¬ 
face  temperatures  with  the  limited  data  available  show  that  the  correct  trends 
and  salient  features  are  predicted.  In  most  instances,  predicted  temperatures 
are  within  a  few  degrees  Celsius  of  measured.  Preliminary  model  performance 
indicates  that  the  technique  incorporates  much  of  the  correct  physics  and  is 
capable  of  being  developed  into  a  reliable,  predictive  model  (Hodge  1989). 

Mathematical  Model 

129.  An  examinati.  n  of  important  physical  phenomena  that  determine  the 
temperature  of  various  bulk  storage  tank  surfaces  is  a  necessary  first 
step  in  devising  an  adeq  ate  predictive  simulation.  Figure  11,  a  cross  sec¬ 
tion  of  a  POL  storage  tank  and  its  surroundings ,  schematically  illustrates  the 
pertinent  energy  transfer  mechanisms.  The  tank  is  partially  filled  with  some 
liquid,  typically  a  petroleum-based  product.  An  air  space  exists  between  the 
liquid  and  the  tank  top.  The  tank  surface  temperatures  are  determined  by  the 
response  of  the  tank,  liquid,  and  air  masses  to  the  energy  transferred  to  or 


\ DIRECT  SOLAR 
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CONVECTION  FROM  TOP 


Figure  11.  POL  storage  tank  energy  transfer  mechanisms 
(SW  -  short-wave,  LW  -  longwave) 
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from  the  masses.  Energy  transfer  within  the  tank's  metal  structure  is  by 
conduction.  Energy  transfer  to  the  tank's  outside  surfaces  is  by  convection 
and  radiation,  and  energy  transfer  to  the  air  and  liquid  contained  within  the 
tank  is  by  convection.  Heat  transfer  to  or  from  the  soil  is  by  conduction. 

130.  Models  of  varying  sophistication  can  be  developed  to  describe  the 
time -dependence  of  the  surface  temperatures  of  the  tank.  Generally,  the  less 
sophisticated  the  model,  the  greater  the  empiricisms  required  and  the  more 
significant  (and  limiting)  the  assumptions  involved.  However,  simple  models 
that  emulate  the  correct  physics  are  capable  of  yielding  quantitatively  accu¬ 
rate  results  (Stoecker  1980) .  An  examination  of  models  of  varying  levels  is 
provided  below. 

131.  Perhaps  the  most  sophisticated  modeling  approach  for  this  problem 
would  be  to  recognize  its  true  3-D,  time -dependent  nature  and  to  model  the 
conduction  through  the  tank  structure  as  well  as  the  ground  using  a  3-D,  time- 
dependent  heat  conduction  computer  code.  In  such  an  effort,  the  3-D,  time- 
dependent  Navier- Stokes  equations,  with  a  turbulence  model  if  needed,  would  be 
used  to  resolve  the  motion  of  the  liquid  and  air  (both  internal  and  external) 
to  define  the  convective  heat  transfer  contributions  to  the  tank  structure  and 
ground.  The  radiative  heat  transfer  contributions  would  similarly  be  defined 
and  resolved.  Such  an  approach,  which  would  require  considerable  computer 
resources,  could  be  expected  to  model  the  process  from  a  refined-as-possible 
theoretical  basis,  but  the  cost  would  be  prohibitive  for  routine  application. 
Such  a  model  would  not  be  easy  to  calibrate  or  to  use  since  model  output  would 
be  much  more  detailed  (i.e.,  spatially  and  temporally)  than  any  feasible  set 
of  experimental  results.  This  approach  is  thus  not  tenable  for  target  model¬ 
ing  requirements. 

132.  An  alternative  to  the  above  approach  is  to  make  extensive  use  of 
lumped- parameter  models  based  on  salient  physics.  Following  such  an  approach, 
four  features  must  be  modeled:  (a)  the  air  and  liquid  masses  inside  the  tank, 
(b)  the  influence  of  the  ground  on  the  liquid  mass,  (c)  the  response  of  the 
tank  structure,  and  (d)  the  convective  and  radiative  heat  transfer  to  the 
surfaces  of  the  tank. 

133.  The  air  and  liquid  masses  within  the  tank  must  be  considered 
separately  because  of  the  great  disparity  in  thermal  capacity  per  unit  volume 
(air  is  order  of  103  J  m*3  “K*1  and  oil  is  order  of  106  J  m'3  °K'X)  .  Thus, 
during  the  typical  diurnal  cycle,  the  liquid  temperature  would  not  vary  more 
than  a  few  degrees,  while  the  inside  air  temperature  would  vary  many  degrees. 
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The  net  result  is  that  the  liquid  tends  to  act  as  a  constant  temperature 
"sink"  for  the  process  while  the  inside  air  temperature  follows  (lags)  the 
tank's  surface  temperatures.  A  reasonable  modeling  procedure  is  to  treat  the 
liquid  and  inside  air  as  single  masses,  each  with  a  single  (average)  tempera¬ 
ture.  The  liquid  and  inside  air  temperatures  would  vary  as  prescribed  by  the 
thermal  capacity  and  the  net  heat  transfer  to  or  from  the  individual  masses. 
For  these  lumped-mass  models,  the  convective  heat  transfer  is  based  on  a  con¬ 
vective  heat  transfer  coefficient  with  the  driving  potential  being  the 
difference  between  the  lumped-mass  temperature  and  the  surface  temperature. 
Existing  or  modified  convective  heat  transfer  correlations  must  be  provided, 
as  the  lumped-mass  approach  does  not  define  the  heat  transfer  (as  does  the 
Navier-Stokes  equations  approach). 

134.  The  ground  provides  additional  mass  that  is  available  for  energy 
interchange  with  the  liquid  (and  outside  air) .  Consistent  with  the  lumped- 
mass  assumption  made  for  the  air  and  liquid  within  the  tank  is  the  assumption 
that  the  energy  transfer  from  the  bottom  of  the  tank  to  the  soil  is  one- 
dimensional  (i.e.,  constant  across  the  bottom  of  the  tank)  and  that  the  con¬ 
duction  within  the  soil  is  also  one -dimensional.  Effects  of  the  soil  heat 
conduction  other  than  directly  below  the  tank  are  ignored. 

135.  The  major  objective  of  this  model  is  the  prediction  of  the  temper¬ 
ature  on  the  surfaces  of  POL  bulk  storage  tanks.  Because  of  the  high  thermal 
conductivity  (and  consequently  low  thermal  resistance) ,  the  tank  walls  can  be 
treated  in  a  lumped-mass  fashion  similar  to  that  of  the  inside  air  and  liquid. 
However,  because  of  the  low  Biot  number  (Bi  <  0.02)  associated  with  the  tank 
walls,  this  is  a  much  better  assumption  for  the  walls  than  for  the  inside  air 
and  liquid  (Kreith  and  Bohn  1986).  Logic  dictates  that  at  least  the  tank  top, 
tank  bottom,  and  walls  adjacent  to  the  inside  air  and  walls  adjacent  to  the 
liquid  be  considered  separately.  Additionally,  because  of  the  diurnal  motion 
of  the  sun,  some  distinction  as  to  tank  wall  orientation  must  be  made.  The 
simplest  viable  division  is  to  consider  east,  south,  west,  and  north  faces 
adjacent  to  the  inside  air  and  liquid;  thus,  eight  wall  portions  and  the  tank 
top  are  considered.  Each  individual  wall  face  is  considered  as  isolated  from 
the  others  so  that  no  conductive  path  exists  between  tank  wall  faces. 

136.  The  remaining  feature  to  be  considered  involves  the  convective  and 
radiative  heat  transfer  to  the  outside  surface  of  the  tank  (Incropera  and  de 
Witt  1981,  Poulikakos  1985,  Sparrow  and  Prata  1985,  and  Young  and  Ozel  1985). 
As  with  the  inside  air  and  liquid,  the  convective  heat  transfer  can  be 
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specified  by  existing  heat  transfer  correlations,  with  the  driving  potentials 
being  the  differences  between  the  outside  ambient  temperature  and  the  surface 
temperatures.  Because  of  the  possibility  of  horizontal  winds  as  well  as 
vertical  buoyancy  forces,  mixed  free  and  forced  convection  effects  must  be 
accounted  for  on  the  exterior  surfaces.  For  this  problem,  the  dominant  energy 
transfer  mechanism  is  the  solar-driven  radiative  heat  transfer.  Shortwave  and 
longwave,  direct  and  diffuse,  contributions  as  well  as  thermal  radiation  from 
the  ground  and  tank  surfaces  must  be  considered.  Because  of  the  sun  path, 
various  surfaces  will  be  exposed  to  varying- intensity  solar  radiation  as  the 
day  progresses  (Weiss  and  Scoggins  1989) . 

137.  Energy  balances  must  be  written  for  each  of  the  12  (tank  top,  tank 
bottom,  eight  wall  segments,  inside  air,  inside  liquid)  lumped-mass  components 
introduced  earlier.  Each  energy  balance  is  based  on  conservation  of  energy  as 
applied  to  a  control  volume;  the  general  form  is  (Hodge  1989) 

-pcVol  =  Eq  (105) 


where 

pcVol  -  thermal  capacity  of  mass  within  control  volume,  W-sec 
T  -  mass  temperature,  °C 
t  -  time,  sec 

Eq  -  summation  of  energy  fluxes  to  and  from  control  volume,  W 
The  vertical  wall  target  surface  temperature  model  (Weiss  and  Scoggins  1989) 
was  used  to  calculate  the  radiative  flux  on  the  fluid- filled  target  walls  and 
to  calculate  the  resulting  physical  temperature  and  radiance  of  the  target 
surfaces.  In  this  way,  the  radiance  of  POL  tanks  is  predicted  in  terms  of 
local  weather  conditions,  terrain  surface  type,  target  surface  emissivity  and 
absorptivity,  latitude,  and  time  of  day. 
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PART  VI:  RADAR  SIGNATURES  OF  TERRAIN  BACKGROUNDS 


138.  Knowledge  of  the  interaction  of  electromagnetic  waves  with  terrain 
features  is  important  to  discriminate  targets  from  background  features  using 
active  or  passive  radar  sensors.  An  advantage  of  using  millimeter  or  micro- 
wave  radar  systems  as  target  locators  is  that  they  function  reasonably  well 
under  adverse  weather  conditions  such  as  cloud  cover,  mist,  or  fog  where  IR 
and  visible  light  sensors  lose  their  effectiveness.  The  main  difficulty  with 
using  radar  detectors  is  the  discrimination  of  target  signals  from  the  back¬ 
ground  clutter;  often  the  background  clutter  masks  the  target  signal  or  causes 
false  target  alarms.  Therefore,  it  is  necessary  to  be  able  to  predict  the 
magnitude  of  the  background  clutter  for  various  types  of  terrain  textures  such 
as  soils,  plowed  fields,  grass,  forest  canopies,  and  urban  areas.  The  rela¬ 
tive  magnitude  of  a  target  signal  and  the  background  clutter  depends  on  the 
operational  characteristics  of  a  radar  transceiver  system,  such  as  the  fre¬ 
quency  and  polarization  of  the  transmitted  and  received  waves,  and  on  the 
depression  angle  of  the  transmitted  waves. 

139.  The  radar  return  signal  from  real  terrain  is  generally  due  to  the 
processes  of  reflection  and  scattering.  Reflection  is  a  surface  effect  while 
scattering  occurs  at  rough  surfaces  or  within  the  volume  of  the  terrain  tex¬ 
ture  material.  The  volume  scattering  is  generally  associated  with  the  random 
variation  of  the  structure  of  terrain  materials  about  some  average  value. 

Both  surface  and  volume  scattering  are  important  for  natural  terrain  materials 
(Ftmg  and  Ulaby  1983)  and  can  be  described  by  differential  scattering  cross 
sections  designated  by  aB  and  ov  ,  respectively.  The  total  differential 
scattering  cross  section  is  written  as 


o  =  at  +  ov 


(106) 


and  the  average  power  received  by  an  antenna  of  an  active  radar  sensor  system 
is  written  as  (Fung  and  Ulaby  1983) 


Px  *  XVUk)3  PtGeGr  fo/R 4  dS 


(107) 
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where 


Pr  -  received  power,  W 

A  -  wavelength  of  electromagnetic  radiation,  m 

Pt  -  transmitted  power,  W 

Gt  -  gain  of  transmitting  antenna 

Gc  —  gain  of  receiving  antenna 
a  -  differential  cross  section,  m2  m'2 

R  -  distance  between  radar  antenna  and  the  scattering  point  on  the 
terrain  surface ,  m 

S  -  area  of  radar  footprint  on  terrain  surface,  m2 

140.  From  Equation  107  it  is  clear  that  the  terrain  characteristics 
enter  the  received  power  calculation  through  the  values  of  the  differential 
cross  section.  The  differential  cross  section  for  the  scattering  of  electro¬ 
magnetic  waves  from  terrain  elements  depends  on  a  number  of  physical 
parameters : 

&.  Surface  microroughness  -  small-scale  terrain  features. 

b.  Surface  mesoroughness  -  medium-scale  terrain  features. 

£,.  Wavelength  of  electromagnetic  radiation. 

d.  Depression  angles  of  receiver  and  transmitter. 

e.  Polarization  of  the  transmitted  and  received  signals- -W,  VH, 
HV,  HH. 

£.  Electrical  properties  of  terrain  elements- -average  values, 
standard  deviations,  and  correlation  lengths  for  the  dielec¬ 
tric  constant,  magnetic  permeability,  and  electrical 
conductivity. 

The  electrical  properties  of  terrain  elements  are  determined  by  soil  and  vege¬ 
tation  type,  moisture  content  of  soil  and  vegetation  type,  moisture  content  of 
soil  and  vegetation,  wavelength  of  the  electromagnetic  waves,  and  soil  and 
vegetation  temperature. 

141.  Many  models  of  varying  complexity  are  available  that  can  represent 
the  differential  scattering  cross  section  of  terrain  textures  in  terms  of  the 
physical  parameters  listed  in  the  above  paragraph  (Fung  and  Ulaby  1983) .  The 
model  used  in  this  report  to  create  the  Regional  Radar  Scattering  Map  is  a 
random  medium  volume  and  rough  surface  radar  wave  scattering  model  for  vege¬ 
tation,  soils,  and  snow  cover  (Zuniga,  Habashy,  and  Kong  1979;  Lee  and  Kong 
1985a, b;  Borgeaud,  Shin,  and  Kong  1987). 

142.  A  three -layer  random  medium  model  is  developed  to  study  the  fully 
polar ime trie  scattering  properties  of  earth  terrain.  The  top  layer  is  modeled 
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as  an  isotropic  random  medium,  the  middle  layer  as  an  anisotropic  random 
medium,  and  the  bottom  layer  as  a  homogeneous  half-space.  Volume  scattering 
effects  of  both  random  media  are  characterized  by  correlation  functions  in 
which  variances  and  correlation  lengths  describe  strengths  of  permittivity 
fluctuations  and  physical  sizes  of  embedded  inhomogeneities,  respectively. 

The  anisotropic  effect  of  the  middle  layer  is  attributed  to  specific  structure 
and  alignment  of  the  scatterers.  With  the  strong  fluctuation  theory,  the  mean 
fields  in  the  random  media  are  derived  under  the  bilocal  approximation  with 
singularities  of  the  dyadic  Green's  functions  properly  taken  into  consider¬ 
ation.  With  the  discrete  scatterer  concept,  effective  permittivities  of  the 
random  media  are  calculated  by  two-phase  mixing  formulas.  Then,  the  distorted 
Born  approximation  is  used  to  calculate  the  covariance  matrix,  which  describes 
the  fully  polarimetric  scattering  properties  of  the  terrain  and  is  used  in 
radar  image  simulation  and  earth  terrain  identification  and  classification 
(Lim  et  al.  1989,  Kong  et  al.  1990). 

143.  A  two- layer  random  medium  model  has  been  successfully  applied  to 
polarimetric  remote  sensing  of  earth  terrain  such  as  vegetation,  meadow,  and 
ice  layer.  The  results  obtained  with  the  three -layer  configuration  have  the 
capability  of  accounting  for  polarimetric  scattering  from  earth  terrain  under 
the  effects  of  weather,  seasonal  variation,  and  atmospheric  conditions  such  as 
forest  under  mist,  meadow  under  fog,  and  ice  under  snow.  The  effects  on 
polarimetric  wave  scattering  due  to  the  top  layer  are  identified  by  comparing 
the  three -layer  model  results  with  those  obtained  from  the  two- layer  model. 

The  enhancement  of  the  radar  returns  due  to  dry- snow  cover  on  top  of  first- 
year  sea  ice  observed  in  the  experimental  data  can  be  explained  using  the 
three -layer  random  medium  model.  The  theoretical  results  are  illustrated  by 
comparing  the  calculated  covariance  matrices  with  the  polarimetric  measurement 
data  (Zuniga,  Habashy,  and  Kong  1979). 

144.  The  scattering  of  electromagnetic  waves  from  a  randomly  perturbed 
periodic  surface  is  formulated  by  the  extended  boundary  condition  (EBC)  method 
and  solved  by  the  small  perturbation  method  (SPM).  The  scattering  from  peri¬ 
odic  surface  is  solved  exactly,  and  this  solution  is  applied  in  the  SPM  to 
solve  for  the  surface  currents  and  scattered  fields  up  to  the  second  order. 

The  random  perturbation  is  modeled  as  a  Gaussian  random  process.  The  theoret¬ 
ical  results  are  illustrated  by  calculating  the  bistatic  and  backscattering 
coefficients.  It  is  shown  that  as  the  correlation  length  of  the  random  rough¬ 
ness  increases,  the  bistatic  scattering  pattern  of  the  scattered  fields  shows 
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several  beams  associated  with  each  Bragg  diffraction  direction  of  the  periodic 
surface.  When  the  correlation  length  becomes  smaller,  the  shape  of  the  beams 
becomes  broader.  The  results  obtained  using  the  EBC/SFtf  method  also  compare 
with  the  results  obtained  from  the  Kirchhoff  approximation.  The  Kirchhoff 
approximation  results  agree  quite  well  with  EBC/SPM  method  results  for  the  HH 
and  W  polarized  backscattering  coefficients  for  small  angles  of  incidence. 
However,  the  Kirchhoff  approximation  does  not  give  depolarized  returns  in  the 
backscattering  direction,  whereas  the  results  obtained  using  the  EBC/SPH 
method  give  significant  depolarized  returns  when  the  incident  direction  is  not 
perpendicular  to  the  row  direction  of  the  periodic  surface  (Yueh,  Shin,  and 
Kong  1988) . 

145.  We  have  studied  the  Mueller  matrix  and  polarization  covariance 
matrix  for  polarimetric  radar  systems.  The  clutter  is  modeled  by  a  layer  of 
random  permittivity,  described  by  a  3-D  correlation  function,  with  variance, 
and  horizontal  and  vertical  correlation  lengths.  This  model  is  applied,  using 
the  wave  theory  with  Born  approximations  carried  to  the  second  order,  to  find 
the  backscattering  elements  of  the  polarimetric  matrices.  It  is  found  that 

8  of  16  elements  of  the  Mueller  matrix  are  identically  zero,  corresponding  to 
a  covariance  matrix  with  four  zero  elements.  Theoretical  predictions  are 
matched  with  experimental  data  for  vegetation  fields  (Borgeaud,  Shin,  and  Kong 
1987). 

146.  The  vegetation  canopy  and  snow-covered  ice  field  have  been  studied 
with  a  three- layer  model,  an  isotropic  random  medium  layer  overlying  an  aniso¬ 
tropic  random  medium.  We  have  calculated  the  dyadic  Green's  functions  of  the 
three- layer  medium  and  the  scattered  electromagnetic  intensities  with  the  Born 
approximation.  The  backscattering  cross  sections  are  evaluated  for  active 
microwave  remote  sensing.  The  theoretical  approach  can  be  extended  to  derive 
the  bistatic  scattering  coefficients  (Borgeaud,  Shin,  and  Kong  1987). 
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PART  VII:  DEVELOPMENT  OF  LASER  SCATTERING  MODELS 


Background 

147.  A  parallel  effort  to  develop  laser  scattering  models  to  describe 
and  predict  the  interaction  of  laser  light  with  natural  terrain  backgrounds 
was  initiated  in  January  1987.  Development  of  a  1-D  semi-empirical  laser 
scattering  model  was  initiated  at  WES  by  the  members  of  the  Background  Signa¬ 
tures  Team.  A  more  comprehensive  and  ambitious  theoretical  model  was  con¬ 
tracted  for  development  (Ishimaru  1978,  Ma  1989). 


Purpose  and  Scope 


148.  The  objective  of  this  research  was  to  develop  comprehensive  and 
practical  theoretical  and  semi -empirical  models  and  computer  codes  for  laser 
scene  generation.  Achieving  this  objective  meant  conducting  extensive 
research  in  the  following  areas: 

a.  Development  of  a  systematically  comprehensive  laser  scattering 
model  based  on  experience  and  research  in  the  theories  of  wave 
scattering  in  random  media.  These  scattering  models  are 
designed  to  predict  the  bidirectional  reflectance  for  terrain 
textures  such  as  soil,  rocks,  grass,  trees,  urban  areas,  agri¬ 
cultural  areas,  and  other  commonly  occurring  background  types. 
These  scattering  models  should  also  predict  the  bidirectional 
reflectance  in  terms  of  frequency,  polarization,  coherence, 
view  angle,  source  angle,  water  content,  weather,  and  season. 
These  models  are  to  be  incorporated  into  a  coarse -resolution 
regional  map  display  and  a  fine -resolution  laser  scene 
display. 

b.  Development  of  a  synthetic  laser  scene  generation  capability 
that  will  predict  what  a  terrain  background  will  look  like  to 
a  laser  scanner.  This  includes  predictions  of  the  correlation 
length  and  standard  deviation  of  the  spatial  fluctuations  of 
the  scattered  intensity  for  each  type  of  terrain  background 
texture . 

c.  Development  of  theoretical  models  to  estimate  the  statistical 
variations,  spatial  correlations,  and  standard  deviations  of 
the  scattered  radiance  from  terrain  in  order  to  aid  detection 
and  identification  of  targets  for  the  automatic  target  detec¬ 
tion  device  under  various  weather  conditions. 

Incorporation  of  useful  theoretical  models  to  account  for  the 
atmospheric  blurring  effects  due  to  dust,  smoke,  aerosols, 
clouds,  and  snow,  including  the  modulation  transfer  function 
and  the  effects  on  spatial  resolution  and  spatial  frequencies. 
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£.  Generation  of  useful  and  practical  computer  codes  and  graphics 
that  can  be  used  in  the  scene  generation  development  program, 
data  interpretations,  and  scene  simulations,  and  regional  map 
displays  of  characteristic  laser  scattering  intensities  as  a 
function  of  terrain  texture. 

f.  Development  of  data- calibrated,  semi-empirical  models  to 

describe  laser  scattering  from  natural  and  man-made  surfaces 
to  meet  near-term  predictive  requirements. 

Semi -empirical  Laser  Scattering  Model 

149.  The  scattering  properties  of  a  rough  terrain  surface  depend  to  a 
large  degree  on  the  following  factors:  laser  altitude  angle,  sensor  depres¬ 
sion  angle,  average  terrain  roughness  amplitude,  and  average  terrain  roughness 
wavelength  (Bass  and  Fuks  1979,  Mahrer  1982,  Carroll  1982,  and  Weiss  and  Scog¬ 
gins  1986) .  These  factors  will  determine  the  fraction  of  the  total  terrain 
surface  area  that  is  illuminated,  the  fraction  of  the  total  terrain  surface 
area  that  is  viewed  by  the  sensor,  and  the  fraction  of  the  viewed  area  that  is 
illuminated  or  shaded.  The  fraction  of  the  viewed  area  that  is  illuminated 
and  shaded  combined  with  the  physical  properties  of  the  illuminated  and  shaded 
areas  is  used  to  calculate  the  scattering  properties  of  rough  terrain  (Weiss 
and  Scoggins  1986). 

150.  Initially,  the  terrain  surface  is  assumed  to  have  a  sinusoidal 
variation  of  elevation.  A  first  step  toward  calculating  the  scattering  prop¬ 
erties  of  rough  terrain  is  to  determine  the  fraction  of  a  sinusoidal  surface 
that  is  illuminated.  Construction  of  a  sensor  view  model  required  the  devel¬ 
opment  of  two  capabilities: 

a.  The  prediction  of  the  fraction  of  the  rough  terrain  area  that 
is  viewed  by  a  sensor  for  a  specified  depression  angle  (view 
altitude  angle) . 

fe.  The  prediction  of  the  fraction  of  the  viewed  area  that  is 
illuminated. 

The  fraction  of  the  viewed  area  that  is  lighted  and  shaded  has  the  greatest 
influence  in  determining  the  laser  scattering  properties  for  rough  terrain. 

151.  An  analysis  similar  to  one  done  to  calculate  the  percentage  of 
illuminated  terrain  area  can  be  used  to  calculate  the  fraction  of  a  rough 
terrain  surface  viewed  by  a  sensor.  The  important  physical  quantity  to  calcu¬ 
late  for  evaluating  sensor  performance  is  the  fraction  of  the  viewed  area  that 
is  illuminated.  This  fraction  determines  the  effective  scattering  properties 
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of  rough  terrain  that  a  sensor  will  detect  in  its  field  of  view  and  the  bi¬ 
directional  reflectance  for  the  case  of  a  laser  source  and  detector. 

152.  Rough  terrain  surfaces  exhibit  little  specular  reflection;  there¬ 
fore,  a  more  realistic  scattering  formalism  must  be  developed.  The  light¬ 
scattering  model  described  in  this  paper  can  be  applied  to  laser  scattering 
from  rough  surfaces  if  the  wavelength  of  the  laser  light  is  much  smaller  than 
the  wavelength  of  the  terrain  elevation  variation  and  also  if  the  terrain 
roughness  wavelength  is  less  than  the  contact  area  of  the  laser  beam  with  the 
terrain  surface.  In  this  case,  the  complex  scattering  process  can  be 
described  by  the  bidirectional  reflectance  function.  For  a  surface  with  no 
hidden  areas,  it  is  assumed  that  the  scattering  is  purely  Lambertian. 

153.  The  rectilinear  propagation  model  presented  below  describes,  in 
general  terms,  the  interaction  of  light  with  rough  terrain  surfaces  (Weiss  and 
Scoggins  1986;  Meeker,  Scoggins,  and  Weisr  1987;  Meeker  1988).  For  laser 
scattering,  this  model  predicts  the  effective  scattering  properties  of  rough 
terrain  in  terms  of  laser  altitude  angle,  sensor  depression  angle,  and  the 
terrain  roughness  amplitude  and  wavelength. 

154.  Based  on  a  1  sinusoidally  varying  surface,  the  semi -empirical 
descriptive  model  calculates  the  fraction  of  the  surface  illuminated  by  the 
radiative  source  for  a  given  incident  angle.  For  a  given  view  angle,  the 
amount  of  illuminated  area  is  determined,  and  the  percentage  of  reflected 
radiation  detected  is  calculated.  The  set  parameters  used  to  generate  the 
descriptive  curves  can  be  broken  down  into  three  groups: 

a.  The  amount  of  backward  scattering  is  depicted  through  the  use 
of  the  fraction  illuminated  calculation,  and  the  forward  scat¬ 
tering  component  is  related  by  the  fraction  viewed 
calculation. 

fe.  A  term  for  diffuse  reflection  includes  both  the  view  angle  and 
the  incident  angle,  while  the  specular  reflection  term  is 
dependent  only  on  incident  angle. 

£.  The  backscatter  and  forward  scatter  segment  of  the  model  has  a 
set  of  parameters  to  weight  the  scattering  contributions 
appropriately,  as  determined  by  the  experimental  data  for  the 
particular  class  of  material. 

In  addition,  the  three  main  groups  all  possess  weighting  coefficients  of  a 
similar  nature  which  are  related  to  each  other  in  such  a  way  that  their  sum  is 
always  equal  to  1. 

155.  Critical  inputs  to  the  model  include  the  wavelength  of  the  inci¬ 
dent  radiation,  the  peak-to-peak  amplitude  of  the  sinusoid  divided  by  its 
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wavelength,  and  the  land  use  classification  index  number.  These  inputs  com¬ 
bine  to  provide  the  model  with  flexibility  in  regard  to  surface  roughness  and 
spectral/reflectance  characteristics  for  various  types  of  soils,  vegetation, 
and  urban  surfaces  necessary  for  an  accurate  prediction  of  background  signa¬ 
tures  for  the  visible  spectra. 

Experimental  Calibration  of  T^gsr  Models 

156.  For  an  experiment  in  conjunction  with  the  National  Aeronautics  and 
Space  Administration's  (NASA)  Goddard  Space  Flight  Center  Code  623  and  300 
laboratories,  a  goniometer  was  adapted  to  accommodate  a  HeNe  laser,  operating 
at  633  nm,  and  fiber-optic  cable  coupled  to  a  silicon  photodiode  detector.  A 
set  of  sample  holders  conforming  to  the  National  Bureau  of  Standards  (NBS) 
specifications  for  bidirectional  reflectance  distribution  function  (BRDF) 
determination  was  constructed.  The  NBS  employs  powdered  Halon  pressed  to  a 
density  of  1  g  cm-3  as  the  reference  for  measuring  BRDF;  therefore,  a  powder 
press  for  Halon  was  constructed  to  NBS  specifications.  Actual  data  values 
were  acquired  through  the  use  of  a  microcomputer  set  at  a  sample  time  of 
either  60  or  400  msec.  All  400-msec  data  were  acquired  using  a  fiber-optic 
cable  with  lower  transmission  values  than  the  cable  used  in  the  60-msec 
measurements . 

157.  Data  were  acquired  at  10-deg*  intervals  from  -70  through  70  deg  in 
the  principal  plane  of  incidence.  In  one  set  of  samples,  the  angular  incre¬ 
ment  was  reduced  to  2.5  deg,  while  some  of  the  data  runs  made  by  NASA  person¬ 
nel  were  in  15-deg  increments.  Because  of  the  configuration  of  the  laser  and 
the  receiver,  direct  backscatter  contributions  could  not  be  measured;  there¬ 
fore,  measurements  were  made  at  5  deg  on  either  side  of  the  incident  angle. 

158.  The  acquired  data  were  then  used  to  calibrate  the  WES  laser 
scattering  model.  By  setting  constants  in  the  model  to  control  the  relative 
contributions  for  various  types  of  scattering  as  a  function  of  the  surface 
roughness  parameters  of  amplitude  and  wavelength,  the  model  output  was 
adjusted  to  resemble  the  data  curves.  Thus,  a  semi -empirical  method  for  the 
determination  of  the  laser  scattering  characteristics  function  was 
constructed. 


*  To  convert  degrees  (angle)  to  radians,  multiply  by  0.01745329. 
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PART  VIII:  BATT7.EFIELD  AIDS 


159.  A  wide  variety  of  electromagnetic  sensors  now  exist  for  use  in 
locating  targets.  These  sensors  operate  by  discriminating  the  target  electro¬ 
magnetic  signatures  from  the  terrain  background  electromagnetic  signatures. 
Often  the  target  is  camouflaged  so  as  to  blend  into  the  terrain  background. 

The  target  and  terrain  electromagnetic  signatures  may  be  similar  in  a  given 
waveband  but  different  in  other  wavelengths  of  the  electromagnetic  spectrum. 
Very  often,  visual,  infrared,  and  radar  target  locators  will  have  different 
efficiencies  for  detecting  targets  that  are  hidden  in  complex  terrain  back¬ 
grounds.  In  addition,  local  weather  conditions  will  severely  affect  the 
operational  performance  of  some  electromagnetic  sensors.  Therefore,  it  is 
important  to  have  battlefield  aids  tk'*-  will  assist  a  commander  in  the  selec¬ 
tion  of  an  appropriate  electromagnetic  sensor  for  specific  targets,  terrain 
backgrounds ,  and  weather  conditions . 

160.  The  battlefield  aids  must  be  relatively  simple  so  as  to  be  mobile 
and  operational  under  battlefield  conditions.  Therefore,  the  battlefield  aids 
should  be  in  the  form  of  maps  or  images  that  can  be  displayed  on  a  portable 
computer  system;  at  the  same  time,  they  must  be  sufficiently  sophisticated  so 
as  to  utilize  the  characteristics  of  the  target  surface  and  the  local  terrain 
and  weather  conditions  to  make  accurate  multispectral  predictions  of  the 
target/background  signature  contrast.  The  battlefield  aids  described  in  this 
report  are  based  on  the  physical  model  that  describes  the  interaction  of  elec¬ 
tromagnetic  waves  with  a  target,  terrain  background,  and  sensor.  These 
battlefield  aids  are  presented  in  the  form  of  color-coded  maps,  a  synthetic 
scene  generator,  and  a  weather -dependent  camouflage  performance  assessment 
algorithm.  The  battlefield  aids  considered  in  this  report  are  as  follows: 

a.  Regional  elevation,  slope,  soil  type,  and  land  use  map. 

b.  Regional  thermal  radiance  map. 

c.  Regional  radar  cross  section  map. 

d.  Regional  laser  reflectance  map. 

e.  Synthetic  IR  scene  generator. 

f.  Camouflage  Effectiveness  Evaluation  System. 

These  battlefield  aids  allow  the  user  to  estimate  immediately  the  target  loca¬ 
tion  effectiveness  of  a  specific  sensor  for  a  known  type  of  target  located  in 
a  specified  terrain  background. 
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Regional  Signature  Maps 


161.  WES  has  developed  the  capability  to  integrate  thermal  terrain 
signature  models  with  terrain  description  data  so  that  maps  of  signature 
values  may  be  created  for  geographic  regions  of  interest  (Weiss  and  Scoggins 
1989) .  Improvements  have  been  made  to  the  software  to  expand  the  accuracy  of 
the  thermal  map  by  utilization  of  the  WESTHERM  combined  surface  model 
described  earlier.  Also,  signature  models  for  radar  cross  section  and  laser 
reflectance  have  been  incorporated  into  map  creation. 

Background 

162.  Input  databases .  To  produce  a  signature  map  of  a  specific  geo¬ 
graphic  region,  data  describing  the  terrain  surface  shape  and  characteristics 
must  be  supplied.  The  terrain  description  data  common  to  the  signature  maps 
are  terrain  elevation,  land  use,  and  soil  types.  Surface  slope  and  azimuth 
(i.e.,  surface  orientation  relative  to  0  deg  north)  are  also  necessary  and  can 
be  supplied  if  available  or  calculated  from  the  elevation  data.  Thus,  five 
sets  of  regional  description  inputs  are  required.  These  data  are  in  the  form 
of  digital  maps  that  consist  of  2-D  arrays  of  numbers  laid  out  on  a  regular 
grid.  Each  grid  element  is  the  value  of  the  respective  parameter  at  that  geo¬ 
graphical  location.  As  an  example,  elevation  could  be  specified  on  a  50-m 
grid,  which  means  that  the  grid  elements  are  the  terrain  elevation  measured  at 
regular  points  spaced  50  m  apart  horizontally  and  vertically.  Whereas  the 
elevation  grid  represents  actual  values  in  meters  or  feet,  the  land  use  and 
soil  type  grids  will  contain  integer  class  numbers  that  are  systematically 
assigned  to  specific  feature  types,  such  as  land  use  class  4  for  deciduous 
forest.  Figure  12  contains  example  data  for  a  site  in  Fort  Hunter  Liggett, 
California.  Elevation,  land  use,  and  soil  type  for  this  site  were  digitized 
from  maps;  slope  and  azimuth  were  calculated  from  the  elevation  data. 

163.  In  addition  to  the  regional  description  database,  another  database 
of  terrain  physical  properties  must  be  available.  This  database  is  specific 
to  the  particular  type  of  signature  map  desired.  For  thermal  signature  maps, 
soil  thermal  properties  are  one  parameter  that  must  be  specified  for  each  type 
of  soil  in  the  region.  Alternatively,  a  laser  reflectance  map  requires  sur¬ 
face  roughness  for  each  land  use  in  the  land  use  data  set.  These  data  are 
cross-references  to  the  class  numbers  in  the  land  use  and  soil  type  data  files 
so  that,  for  any  location  in  the  region  of  interest,  all  parameters  necessary 
to  execute  the  specific  signature  model  are  known.  Table  1  is  a  summary  of 
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Figure  12.  Geographic  information  system 
database 
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Table  1 

Summary  of  Database  Inputs  for  Signature  Mao  Models 


Thermal 

Thermal  conductivity 
Thermal  diffusivity 

Surface/vegetation  longwave  emissivity 
Surface/vegetation  shortwave  absorptivity 
Surface  moisture  content 
Vegetation  cover  factor 
Vegetation  height 


Radar 

Root -mean- squared  surface  slope 
Complex  dielectric  constant 
Layer  thickness 

Correlation  length  of  dielectric  constant 


Laser 

Characteristic  surface  roughness  dimension 
Surface  roughness  amplitude 

Bidirectional  reflection  distribution  function  (optional) 


the  types  of  input  required  for  the  thermal,  radar,  and  laser  signature 
models . 

164.  One  final  type  of  input  data  is  scenario-specific  data  entered  by 
the  user  at  the  time  of  map  generation.  These  data  also  depend  on  the  partic¬ 
ular  type  of  map  desired.  A  thermal  map  will  require  the  user  to  enter  the 
weather  conditions  near  the  tim^  of  interest.  The  radar  and  laser  maps 
require  sensor  position,  orientation,  and  operational  wavelength. 

165.  Man  generation.  Given  the  presence  of  the  databases  described 
above,  specific  signature  maps  can  be  generated  using  software  developed  at 
WES  (Weiss  and  Scoggins  1989).  Improvements  to  WES  software  have  been  made  to 
better  account  for  surface  slope  and  azimuth  and  to  allow  generation  of  radar 
cross-section  and  laser  reflectance  maps.  The  procedure  by  which  the  software 
produces  the  desired  map  is  described  below. 

&.  Scan  terrain  region  description  database  to  identify  unique 
combinations  of  classes,  called  class  vectors.  This  list  is 
specific  to  the  signature  model  and  constitutes  the  number  and 
type  of  model  executions  for  the  region.  This  is  a  one-time 
operation  given  that  the  region  of  interest  does  not  change. 
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t>.  Repeatedly  execute  the  desired  signature  model  for  the  class 
vectors  found  in  step  a.  These  class  vectors  are  cross- 
referenced  to  the  parameter  database  to  provide  physical  prop¬ 
erties  necessary  for  the  model  to  generate  a  signature  value. 
Create  signature  for  each  class  vector. 

c.  Scan  the  elevation,  land  use,  soil  type,  slope,  and  azimuth 
maps  in  parallel  and  produce  a  class  vector  for  each  grid 
element.  Cross-reference  the  signature  value  to  the  element 
class  vector.  Create  a  new  map  by  building  a  grid  of  signa¬ 
ture  values  that  are  spatially  registered  with  the  terrain 
region  description  files. 

Example  mans 

166.  Figures  13-15  are  examples  of  thermal,  radar  cross-section,  and 
laser  reflectance  maps,  respectively.  These  maps  were  produced  for  the 
Fort  Hunter  Liggett  region  represented  by  the  terrain  description  maps  shown 
in  Figure  12.  The  thermal  map  represents  the  radiance  from  surfaces  for  a 
clear  day  in  September  at  1400  hr.  Weather  inputs  used  to  produce  Figure  13 
are  given  in  Table  2.  The  radar  cross-section  map  (Figure  14)  is  for  the  same 
area  for  a  system  operating  in  the  35 -GHz  frequency  band  oriented  at  45  deg 
relative  to  north  and  with  an  elevation  of  15  deg  from  horizontal.  Finally, 
the  laser  reflectance  map  (Figure  15)  is  for  the  Fort  Hunter  Liggett  area 
illuminated  by  a  laser  operating  at  1.06  n  with  the  same  source  orientation  as 
with  the  radar.  Unlike  the  radar,  the  laser  map  represents  reflectance 
received  by  a  sensor  not  at  the  source  location.  The  sensor  location  for  the 
laser  map  in  Figure  15  is  335  deg  relative  to  north  and  15  deg  relative  to 
horizontal. 


Thermal  Infrared  Computer  Scene  Generation 

167.  Surface  temperature  models  allow  simulation  of  temperatures  for 
major  terrain  surface  types.  The  models  are  implemented  as  subroutines  to  a 
driving  program  to  facilitate  temperature  predictions  of  the  surfaces  on  the 
terrain  map.  A  specific  terrain  map  is  first  scanned  by  a  series  of  prepro¬ 
cessor  programs,  and  a  list  of  terrain  types,  azimuth  angles,  and  slope  angle 
classes  is  written  to  a  file.  Thus,  when  a  scene  is  desired  for  specific 
weather  conditions  and  time  of  day,  this  file  is  read  by  the  model  driver 
program  for  the  necessary  class  combinations. 

168.  Based  on  the  terrain  type  index  (grass,  trees,  etc.),  parameters 
needed  by  the  appropriate  model  are  extracted  from  a  parameter  database  and 
passed  to  the  correct  model  subprogram.  The  resulting  temperatures  are 


70 


Suxfa.ce  Temperature,  Deg  Celsius 
Figure  13.  Regional  thermal  radiance  map 
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Soattsriaf  Cron  Section,  dB 

Figure  14.  Regional  radar  cross-section  map 
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Figure  15.  Regional  laser  reflectance  map 
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Table  2 

Weather  Conditions  Used  for  Regional  Thermal  Map  Generation 


Time 

hr 

Air 

Temperature 

°C 

Relative 

Humidity 

% 

Cloud 

Cover 

% 

Wind 

Speed 

m/sec 

1 

10.6 

89 

0 

0.31 

2 

10.5 

90 

0 

0.67 

3 

9.4 

90 

0 

0.05 

4 

9.5 

91 

0 

0.21 

5 

8.7 

91 

0 

0.62 

6 

8.2 

91 

100 

0.57 

7 

10.8 

91 

100 

0.1 

8 

13.0 

87 

100 

0.26 

9 

15.2 

77 

50 

0.98 

10 

19.1 

55 

0 

1.34 

11 

21.3 

43 

0 

1.59 

12 

23.7 

34 

0 

1.39 

13 

25.2 

21 

0 

1.95 

14 

26.7 

21 

0 

1.49 

15 

7.6 

18 

0 

1.39 

16 

27.4 

17 

0 

2.36 

17 

26.1 

20 

0 

3.03 

18 

24.3 

23 

0 

1.18 

19 

19.6 

41 

0 

0.26 

20 

18.5 

45 

0 

1.13 

21 

16.3 

54 

0 

0.72 

22 

15.2 

59 

0 

0.51 

23 

12.7 

75 

0 

0.36 

24 

12.3 

81 

0 

0.21 

converted  to  blackbody  radiances  by  an  approximation  equation  (Weiss  and 
Scoggins  1989) .  Values  for  the  standard  deviation  and  correlation  length  of 
radiance  are  selected  from  a  list  derived  by  analysis  of  actual  thermal  imag¬ 
ery  of  the  particular  texture  type.  Finally,  the  mean,  standard  deviation, 
and  correlation  length  of  the  radiance  for  each  surface  type  are  written  to  an 
intermediate  file  for  use  by  the  scene  generator  software. 
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169.  Such  a  scene  generator  Is  in  use  at  WES  for  creation  of  simulated 
thermal  IR  scenes  (Weiss  and  Scoggins  1989) .  Modules  that  degrade  the  "per¬ 
fect"  scene  output  by  the  scene  generator  software  have  been  added  to  include 
the  effects  of  the  thermal  sensor  so  that  the  final  scene  will  more  accurately 
reflect  what  would  be  displayed,  such  as  on  a  monitor. 

Scene  Generator  Program 


Background 

170.  The  structures  of  the  objects  and  terrain  backgrounds  in  the 
database  which  serve  as  input  to  the  scene  generator  are  approximated  by  tri¬ 
angular  facets.  The  computer  program  reads  the  coordinates  of  all  three  ver¬ 
tices  of  the  triangular  facets  from  a  file,  along  with  the  direction  cosines 
of  the  normal  of  each  facet.  The  direction  cosines  of  the  facets  will  be  used 
by  the  program  to  determine  which  facets  are  seen  from  the  sensor's  point  of 
view.  Background  surface- type  information  is  derived  from  a  digital  terrain 
map.  The  terrain  database  is  an  encoding  of  information  regarding  the  place¬ 
ment  of  background  texture  types.  Finally,  terrain  surface  boundary  data, 
such  as  road  edges,  are  inserted  into  the  terrain  map.  These  data  are  digi¬ 
tized  from  conventional  maps  and  are  merged  in  such  a  way  that  boundary  coor¬ 
dinate  data  may  be  at  a  higher  resolution  than  terrain  facets.  This  final 
terrain  map  is  the  geometric  model  used  by  the  scene  generator  for  background 
scene  production  (Weiss  and  Scoggins  1989) . 

171.  To  construct  a  scene,  the  detector  position,  viewing  direction, 
and  orientation  are  specified  in  a  data  file.  The  angular  field  of  view  of 
the  detector  determines  the  footprint  or  scene  size  on  the  ground  plane  which 
is  centered  about  the  line  of  sight  between  the  sensor  and  ground  plane.  The 
x,  y,  and  z  coordinates  of  the  sensor  relative  to  the  origin  of  the  terrain 
database  must  be  entered  into  the  scene  generator  program.  In  addition,  the 
aspect  angle,  azimuth  angle,  and  depression  angle  of  the  sensor  must  be  speci¬ 
fied.  These  parameters,  along  with  sensor  resolution  and  the  size  of  the 
desired  scene  in  number  of  pixels  and  rows,  are  entered  into  a  file  which  is 
then  read  by  the  scene  generator.  This  procedure  defines  one  frame  for  the 
scene  generator  to  create.  The  software  is  structured  so  that  the  parameter 
file  is  repeatedly  read  for  additional  frames  for  the  scene  generator  to  cre¬ 
ate  until  an  end  of  file,  at  which  point  the  program  terminates. 
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172.  The  "brightness"  of  each  scene  pixel  is  related  to  the  energy 
emitted,  or  radiance,  from  the  terrain/object  surfaces.  Each  facet  is 
assigned  a  radiance  value  based  on  its  texture  type,  which  has  been  calculated 
by  the  surface  temperature  models,  and  standard  deviations  and  correlation 
lengths  that  describe  the  variability  of  the  texture  type.  In  general,  the 
intrinsic  radiance  of  a  surface  depends  on  the  orientation  relative  to  the 
sun,  the  surface  characteristics  such  as  absorptivity  and  roughness,  and  the 
thermal  inertia  of  the  subsurface  material. 

Sensor  model 

173.  The  sensor  model  degrades  the  synthetic  image  created  by  the  scene 
generator,  accounting  for  the  effects  of  sensor  optics,  detector  response,  and 
analog- to-digital  conversion.  The  process  consists  of  three  main  steps: 
transfer  of  the  radiant  energy  flux  through  the  sensor  optics,  conversion  of 
the  energy  to  voltage  by  the  detector  including  scanning  effects,  and  digiti¬ 
zation  of  the  voltage  signal.  Each  of  these  steps  removes  detail  from  the 
perfect  scene  produced  by  the  scene  generator. 

174.  Modulation  transfer  function.  The  modulation  transfer  function 
accounts  for  the  blurring  effect  of  the  sensor  optics  on  the  perfect  image  and 
calculates  the  irradiance  incident  on  the  detector  element  from  each  input 
pixel  (Hudson  1969) .  Any  optical  system  with  a  finite  size  aperture  will  blur 
a  point  source  so  that  the  image  will  be  a  spot  rather  than  a  point.  Also, 
the  irradiance,  or  energy  incident  on  the  detector,  is  a  function  of  the  sen¬ 
sor  field  of  view  as  defined  by  the  optics  focal  length.  To  account  for  this, 
a  convolution  or  modulation  of  the  input  image  with  a  2-D  Gaussian  point- 
spread  function  is  performed.  The  Gaussian  function  is  defined  by  a  mean 
related  to  the  effective  field  of  view  and  variances  in  the  x  and  y  dimensions 
that  incorporate  blurring  effects.  The  calculation  is  given  by 


I\  (x,y) 


£  I. 


x  G(x-l,y-j) 


(108) 


where 

Ix  -  modified  image 

m,n  -  extent  of  convolution  defined  by  variances 
I0  -  original  perfect  image 
G  -  zero  mean  Gaussian  function 


79 


The  mean  variances,  entered  by  the  user,  are  related  to  specific  characteris¬ 
tics  of  the  optics  such  as  the  aperture,  field  stop,  and  lens  diameter  and  are 
generally  provided  by  the  system  manufacturer. 

175.  Detector  voltage  output.  The  next  step  in  the  sensor  model  is  to 
calculate  the  voltage  output  from  the  sensor’s  detector  element  (Hudson  1969). 
This  incorporates  the  responsivity  of  detector  and  scanning  interlace  in  con¬ 
version  of  the  irradiance  to  voltage.  Responsivity  is  a  wavelength -dependent 
constant  for  a  detector  material  that  relates  the  irradiance,  in  watts  per 
square  meter,  to  output  voltage.  Also,  the  noise  imposed  on  the  image  by  the 
detector  can  be  calculated  from  responsivity  and  detector  size. 

176.  Unlike  a  real  image,  which  is  a  continuous  brightness  function, 
interlace  effects  occur  when  the  scanning  system  does  not  exactly  align  the 
detector  element  on  each  pixel  of  the  digital  image  produced  by  the  scene 
generator.  Also,  the  detector -optics  combination  may  view  an  area  greater 
than  one  scene  pixel  in  size.  To  account  for  this  properly,  the  initial  scene 
must  be  generated  at  least  twice  the  resolution  of  the  final  image,  say  at 
200  by  200  pixels  for  a  final  image  of  100  by  100  pixels.  The  higher  resolu¬ 
tion  image  is  then  resampled  and  averaged  to  a  lower  resolution  to  simulate 
this  effect. 

177.  Digital- to-analog  conversion.  In  most  imaging  sensor  systems,  the 
output  voltage  is  converted  to  a  number  for  display  by  computer  on  a  monitor. 
An  automatic  gain  control  is  normally  used  to  adjust  the  dynamic  range  of  the 
detector  voltage  to  match  the  scale  of  the  analog- to -digital  conversion  sys¬ 
tem,  generally  0  to  255.  This  is  done  to  optimize  the  contrast  on  the  image 
finally  displayed  on  the  monitor.  Such  modeling  simply  amounts  to  scaling  the 
values  output  by  the  detector- scanner  simulation  to  the  desired  digital  range. 
Input  specification  and  execution 

178.  To  simulate  a  specific  sensor  system,  details  of  the  system  must 
be  entered  into  the  sensor  simulator.  The  system  currently  has  data  necessary 
to  model  the  Inframetrics  and  AGA  infrared  scanners,  with  facilities  for 
entering  data  to  simulate  other  systems  if  such  data  are  available.  The  scene 
generator  itself  is  executed  initially  to  produce  the  perfect  image,  requiring 
only  the  terrain  and  object  geometry  and  view  specifications  (Weiss  and 
Scoggins  1989) .  The  sensor  model  utilizes  the  output  of  the  scene  generator 
to  produce  a  degraded  image  that  is  a  closer  approximation  to  what  would  actu¬ 
ally  be  seen  on  the  monitor  of  scanning  imager.  Since  the  sensor  model  is  a 
proto-processor,  it  is  possible  to  use  the  software  with  other  scene 
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generators,  provided  that  the  scene  generator  produces  radiance  values  as 
output . 


ramntifl.ape  Effectiveness  Evaluation  System  (GEES') 


Background 

179.  The  CEES  is  a  set  of  software  that  combines  the  thermal  radiance 
prediction  models  for  targets  and  terrain  background  with  the  Night  Vision 
Laboratory  Static  Performance  Model  for  Thermal  Viewing  Systems  (Ratches  et 
al.  1975)  to  obtain  a  model  that  predicts  target  detection  probability  at  a 
specified  range  or  target  detection  range  for  a  0.5  detection  probability 
(Weiss  and  Scoggins  1989).  The  TSTM,  VEGIE,  and  VCTM  thermal  models  are  used 
to  predict  terrain  background  radiance  for  use  with  IR  sensors.  CEES  allows 
interactive  specification  of  values  for  parameters  such  as  emissivity,  as  well 
as  other  parameters  that  influence  the  surface  temperature .  The  software  is 
designed  to  allow  for  easy  operation  by  the  user  and  for  output  of  the  results 
in  tabular  as  well  as  graphic  form  (Analytics,  Inc.  1985).  Primarily  designed 
as  a  desktop  tool  for  analysis  of  multispectral  camouflage  on  fixed  installa¬ 
tions  and  other  high-value  targets,  CEES  has  an  input  interface  to  aid  the 
user  in  the  evaluation  and  selection  of  model  inputs.  The  following  describes 
the  incorporation  of  a  passive  millimeter  wave  (mmw)  receiver  model  into  CEES. 
Technical  approach 

180.  Surface  signature  models.  CEES  uses  the  WES  thermal  models  to 
estimate  background  and  target  temperatures  for  analysis  purposes.  These 
models  are  executed  as  subroutines  in  the  system  and  provide  output  for  use  by 
the  sensor-specific  portions  of  the  systems.  As  described  above,  TSTM,  VEGIE, 
and  VCTM  take  as  input  meteorological  data  commonly  collected  by  standard 
weather  stations.  These  data  are  used  to  calculate  a  time -dependent  surface 
boundary  condition  for  use  with  a  finite-difference  formulation  of  the  thermal 
conduction  equation.  Since  the  models  predict  physical  temperatures,  their 
output  can  be  used  by  both  the  mmw  and  IR  modules  of  the  system  performance 
calculations . 

181.  Passive  millimeter  wave  model.  The  mmw  model  incorporated  into 
the  CEES  program  models  a  passive  mmw  receiver.  A  passive  receiver  detects 
the  reflected  and  thermal  radiation  in  the  mmw  spectral  bank  of  bodies  and 
materials  that  are  at  temperatures  above  absolute  zero.  The  output  of  the  mmw 
radiometer  is  proportional  to  the  intensity  of  observed  radiation  within  the 
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bandwidth  of  the  receiver.  The  model  determines  effective  temperatures  for 
targets  and  backgrounds  at  various  ranges  based  on  (a)  slant  range  to  target, 
(b)  target  and  background  emissivities ,  (c)  reflectance,  and  (d)  atmospheric 
losses  (clear,  fog,  rain). 

182.  The  acquisition  range  is  then  calculated  based  on  a  minimum  dis¬ 
cernible  temperature  difference  between  the  target  and  background.  The  mmw 
model  was  designed  for  an  operating  frequency  of  36  GHz.  However,  it  could 
easily  be  adopted  to  other  frequencies  if  the  emissivity  and  atmospherics 
tables  for  different  materials  at  these  frequencies  were  added. 

183.  The  mmw  model  uses  as  inputs  the  physical  temperatures  for  the 
target  and  background,  calculated  as  described  above  with  the  WES  thermal 
models .  These  are  then  attenuated  through  the  atmosphere  to  ob*  »in  the  effec 
tive  temperatures  at  the  sensor.  The  following  equation  is.  used  to  calculate 
effective  temperatures  while  allowing  for  loss  and  self-emissions  from  the 
atmosphere . 
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where 
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T. 
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effective  temperature 
slant  range 

emissivity  of  target  or  background 
physical  temperature  of  target  or  background 
sky  temperature 
atmospheric  loss  factor 
ambient  temperature 


The  term  is  the  Beer's  Law  term  to  account  for  atmospheric  effects 

in  the  millimeter  wave  spectral  region  (Poradish  1982). 

184.  The  mmw  model  calculates  acquisition  ranges  by  initially  setting 
the  range  to  the  sensor  at  1  km  and  then  increasing  this  distance  in  5-km 
intervals  until  the  target  can  no  longer  be  acquired.  Target  acquisition  is 
determined  by  comparing  the  difference  between  effective  temperatures  of  the 
target  and  background  at  a  given  range  with  the  minimum  discernible  AT  value 
as  input  by  the  user.  This  value  normally  ranges  from  0.1  to  10  °K.  Once  the 
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target  acquisition  range  has  been  determined  to  within  a  5-km  interval,  the 
interval  size  is  reduced  to  1  km,  and  the  search  continues  until  an  acquisi¬ 
tion  range  accurate  to  1  km  1 s  determined. 
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PART  IX:  SUMMARY  AND  RECOMMENDATIONS 


Summary 

185.  Several  multispectral  terrain  background  signature  models  have 
been  described  in  this  report.  These  models  are  based  on  physical  principles 
and  describe  the  interaction  of  electromagnetic  waves  with  terrain  back¬ 
grounds.  The  thermal  models  predict  the  physical  temperature  and  the  radiance 
of  rough  terrain  backgrounds  under  ambient  weather  conditions  for  the  relevant 
3-  to  5-/im  and  8-  to  14-^m  wavebands.  The  radar  scattering  models  calculate 
the  differential  scattering  cross  section  of  a  rough  terrain  surface  in  terms 
of  surface  and  volume  scattering  parameters,  frequency,  polarization  (W,  VH, 
HV,  or  HH) ,  and  depression  angles  of  the  transmitter  and  receiver.  In  a 
similar  manner,  a  laser  scattering  model  is  available  that  predicts  bidirec¬ 
tional  reflectance  and  surface  depolarization  for  the  scattering  of  laser 
light  from  rough  terrain.  These  models  have  been  incorporated  into  battle¬ 
field  aids  useful  in  assisting  battlefield  commanders  to  select  the  appro¬ 
priate  electromagnetic  sensor  system  for  given  target,  terrain,  and  weather 
conditions.  The  models  that  were  developed  are  summarized  below. 

Thermal  signature  models 

186.  Terrain  Surface  Temperature  Model  (TSTM) .  This  model  predicts 
physical  temperature  and  radiance  of  smooth  nonvegetative  surfaces  in  terms  of 
time  of  day,  ambient  weather  conditions,  and  radiative  surface  conditions  such 
as  emissivity  and  absorptivity,  and  in  terms  of  heat  conductivity  and  thermal 
diffusivity  of  the  semi-infinite  half-space  used  to  represent  the  ground.  An 
energy  budget  condition  at  the  ground  surface  is  used  to  determine  the  physi¬ 
cal  temperature  at  the  ground  surface,  and  a  heat  conduction  equation  is 
solved  to  obtain  the  subsurface  temperatures.  The  basic  physical  processes  of 
insolation,  IR  surface  emission,  absorption  of  IR  radiation  from  the  sky,  heat 
conduction  in  the  ground,  convection  heat  loss  to  the  atmosphere,  and  heat 
exchange  due  to  evaporation  and  the  formation  of  dew  (phase  changes)  are  con¬ 
sidered  in  the  surface  energy  budget  condition. 

187.  Two-Dimensional  Terrain  Surface  Temperature  Model  (2D-TSTM).  This 
model  is  a  2-D  version  of  the  TSTM  that  predicts  the  physical  temperature  and 
radiance  of  smooth  nonvegetative  terrain  surfaces  but  also  allows  for  the 
possibility  of  lateral  variation  of  the  thermal  properties  of  the  ground.  The 
physical  principles  for  this  model  are  essentially  the  same  as  for  the  TSTM. 
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188.  Vegetative  Surface  Temperature  Model  (VEGIE) .  This  model  predicts 
the  physical  temperature  and  radiance  of  grass.  It  is  a  subprogram  of  the 
TSTM  and  extends  this  model  by  adding  a  thin  layer  of  vegetation  over  the 
ground  surface.  Thermal  IR  energy  is  reflected  back  and  forth  between  the  two 
surfaces  of  this  layer,  and  energy  budget  conditions  are  imposed  simultane¬ 
ously  at  the  top  and  bottom  surfaces  of  the  layer.  In  this  way,  the  physical 
temperature  and  radiance  at  the  top  surface  can  be  determined  using  the  basic 
principles  of  the  TSTM. 

189.  Vegetation  Canopy  Temperature  Model  (VCTM1 .  This  model  predicts 
the  physical  temperature  and  radiance  of  tree  canopies.  The  model  is  a  plane- 
parallel  abstraction  of  a  vegetation  canopy  that  is  divided  into  three  hori¬ 
zontal  layers.  The  radiation  transfer  equation  is  solved  for  this  system  sub¬ 
ject  to  energy  budget  conditions  at  each  of  the  boundaries.  No  heat  storage 
is  allowed  within  the  canopy.  Detailed  canopy  geometry  information  is 
required  for  each  layer,  including  leaf  area  index  and  leaf  slope  distribu¬ 
tion.  In  addition,  the  emission  and  absorption  characteristics  of  the  ground 
and  canopy  vegetation  and  the  canopy  stomatal  resistance  to  water  vapor  diffu¬ 
sion  are  required  to  be  known  parameters. 

190.  Rough  Terrain  Surface  Temperature  Model  (RTSTM) ,  This  model  is  an 
extension  of  the  basic  TSTM/VEGIE  model  to  the  case  of  a  rough  surface 
described  by  an  average  amplitude  and  a  correlation  length.  The  thermal  radi¬ 
ance  is  calculated  by  evaluating  the  fraction  of  the  area  viewed  by  an  IR 
sensor  that  is  sunlit.  Because  this  fraction  depends  on  both  the  depression 
angle  of  the  sensor  and  the  altitude  angle  of  the  sun,  it  follows  that  the 
predicted  radiance  values  also  depend  on  these  two  angles.  The  standard  input 
parameters  for  the  TSTM/VEGIE  model  are  also  required  to  run  the  RTSTM. 

191.  Vertical  Wall  Target  Surface  Temperature  Model  (VWTSTM) .  This 
model  is  an  extension  of  the  TSTM  to  the  case  of  a  vertical  wall  of  a  solid 
body.  For  the  case  of  a  vertical  wall,  additional  radiative  fluxes  arise  due 
to  reflected  visible  and  IR  radiation  from  the  ground  to  the  wall.  In  addi¬ 
tion,  there  is  the  IR  ground  emission  term  that  impinges  on  the  wall.  The 
model  determines  the  physical  temperature  and  radiance  of  the  vertical  walls 
of  solid  body  targets. 

192.  Fluid-Filled  Target  Surface  Temperature  Model  (FFTSTM).  This 
model  predicts  the  physical  temperature  and  radiance  of  the  vertical  wall  of  a 
fluid- filled  target  such  as  a  POL  tank.  The  FFTSTM  is  essentially  the  VWTSTM 
with  a  convective  heat  transfer  process  replacing  the  heat  conduction 
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mechanism.  The  FFTSTM  can  be  used  to  evaluate  camouflage  applied  to  buildings 
such  as  aircraft  shelters. 

Radar  model 

193.  The  model  predicts  the  differential  scattering  cross  section  of 
vegetative  and  nonvegetative  surfaces  in  terms  of  surface  roughness  parameters 
and  in  terms  of  the  physical  structure  of  the  terrain  elements.  This  model  is 
a  full  wave  solution  of  Maxwell's  equations  for  electromagnetic  waves  scatter¬ 
ing  in  random  media  such  as  soils  and  vegetation.  The  electromagnetic  struc¬ 
ture  of  terrain  materials  is  represented  by  mean  values  of  the  dielectric 
constant  and  the  magnetic  permeability  and  by  the  standard  deviations  and 
correlation  lengths  for  these  quantities.  The  differential  scattering  cross 
section  is  predicted  for  a  frequency  range  of  0  -»  100  GHz  for  any  value  of  the 
depression  angles  of  the  receiver  and  transmitter  and  for  the  four  possible 
transmitter  and  receiver  polarization  modes- -HH,  HV,  VH,  and  W. 

Laser  models 

194.  The  development  of  models  that  describe  the  scattering  of  laser 
light  from  real  terrain  features  such  as  soils,  grass,  leaves,  bushes,  and 
trees  has  proceeded  in  two  directions:  (a)  an  in-house  research  program  to 
develop  a  semi -empirical  laser  scattering  model  and  (b)  a  contracted  research 
program  to  develop  a  rigorous,  physically  based  laser  scattering  model.  The 
semi -empirical  model  is  based  on  the  amount  of  scattering  area  available  in 
the  field  of  view  of  a  laser  sensor.  This  is  calculated  using  a  sinusoidal 
representation  of  the  rough  terrain  surface  and  calculating  the  fraction  of 
the  viewed  area  illuminated  by  the  laser.  In  this  way,  the  bidirectional 
reflectance  is  predicted  in  terms  of  the  depression  angles  of  the  laser  and 
the  detector.  The  more  fundamental  laser  scattering  model  considers  both 
volume  and  surface  scattering  and  assumes  the  soil  and  vegetation  to  be  random 
media.  A  radiative  transfer  assumption  approximates  the  solution  of  Maxwell's 
equations  for  the  scattering  of  light  in  the  media. 

Battlef^ej^ald.s 

195.  The  models  described  earlier  in  this  report  have  been  incorporated 
into  various  visual  displays  that  can  be  used  as  battlefield  aids  to  field 
commanders  who  must  select  the  most  appropriate  weapon  sensor  system  for  known 
target,  background,  and  weather  conditions.  The  battlefield  aids  are  devel¬ 
oped  to  operate  on  PC- type  computers  in  order  to  have  the  mobility  that  field 
units  require. 
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196.  The  battlefield  aids  developed  at  VES  include: 

a.  Regional  elevation,  slope,  and  land  use  maos.  These  maps 
display  the  elevation,  slope  and  type  of  terrain  for  each 
region  within  the  boundaries  of  the  maps  by  using  a  color- 
coded  scheme  to  indicate  elevation,  slope,  and  land  use.  Gen¬ 
erally,  a  100-m  grid  system  is  used,  but  for  some  areas  data 
are  available  on  a  20-m  grid. 

b.  Regional  thermal  radiance  map.  The  terrain  temperature  and 
radiance  prediction  models  TSTM,  VEGIE,  and  VCTM  are  applied 
to  each  area  of  known  elevation,  slope  ,  and  land  use  type,  and 
the  radiance  and  physical  temperature  are  calculated.  A  color 
code  displays  the  thermal  properties  of  each  area  within  the 
map.  These  maps  can  be  generated  for  any  time  of  day  and 
night  and  for  the  weather  conditions  expected  during  the  time 
of  engagement  with  a  target.  This  battlefield  aid  can  select 
the  proper  thermal  camouflage  for  military  equipment  and 
structures.  It  can  also  estimate  the  expected  thermal  con¬ 
trast  an  enemy  target  will  have  at  the  time  of  engagement. 

£.  Regional  radar  cross-section  mao.  The  empirical  terrain  radar 
cross-section  prediction  model  developed  at  WES  has  been  com¬ 
bined  with  the  terrain  land  use  map  to  create  a  color-coded 
radar  cross-section  map  that  displays  the  magnitude  of  the 
differential  scattering  cross  section  of  the  various  terrain 
elements  of  the  background  of  a  target.  In  this  way  the 
expected  conspicuity  of  an  enemy  target  can  be  estimated  for  a 
given  radar  system.  On  the  other  hand,  the  radar  cross- 
section  color  map  can  be  used  to  select  proper  camouflage 
against  aircraft  radars. 

d.  Regional  laser  reflectance  map.  The  WES  semi -empirical  laser 
reflectance  prediction  model  was  joined  to  the  terrain  land 
use  map  in  order  to  produce  a  laser  reflectance  map  that  uses 
a  color-coded  scheme  to  represent  the  magnitude  of  the 
reflectance  for  each  terrain  type  in  the  area  of  interest. 

e.  Synthetic  IR  scene  generation.  The  terrain  elevation,  slope, 
and  land  use  database,  together  with  the  physically  based 
terrain  surface  temperature  prediction  models,  have  been  com¬ 
bined  with  a  scene  generation  computer  program  to  produce  a 
synthetic  IR  target  in  background  scene  generation  capability. 
The  scene  generation  program  consists  of  a  target  geometry 
model  and  a  background  geometry  model  that  construct  the  tar¬ 
get  and  background  In  terms  of  planar  triangular  facets.  The 
angular  field  of  view  of  the  IR  detector  determines  a  foot¬ 
print  on  the  ground  plane  characterized  by  the  triangular 
facets.  Each  triangular  terrain  facet  is  associated  with  a 
terrain  texture  type  and  slope,  and  the  TSTM,  VEGIE,  and  VCTM 
models  are  invoked  to  calculate  the  radiance  of  each  facet. 

The  TSTM  is  used  to  calculate  the  radiance  of  the  target 
facets.  The  radiance  of  each  facet  depends  on  its  orientation 
to  the  sun,  consequently  introducing  a  time  dependence  to  the 
radiance  of  the  facets  as  the  sun  moves  across  the  sky.  Addi¬ 
tional  time  dependence  is  introduced  by  a  varying  air  tempera¬ 
ture.  The  effects  of  weather  are  introduced  through  the  basic 
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TSTM  program,  which  includes  effects  of  cloud  cover  and  cloud 
type.  An  EOSAEL  subprogram  accounts  for  the  attenuation  of 
the  1R  image  as  it  is  projected  through  the  atmosphere.  In 
this  way,  a  time -dependent  and  weather -dependent  synthetic  1R 
image  of  a  target  in  a  realistic  background  can  be  created. 

f.  Camoufl age  Effectiveness  Evaluation  System.  The  CEES  is  a  set 
of  software  that  combines  the  thermal  radiance  prediction 
models  for  targets  and  terrain  backgrounds  with  the  Night 
Vision  Laboratory  Static  Performance  Model  for  Thermal  Viewing 
Systems  to  obtain  a  model  that  predicts  target  detection  prob¬ 
ability  at  a  specified  range  or  target  detection  range  for  a 
0.5  target  detection  probability.  The  TSTM  calculates  the 
target  radiance,  while  the  TSTM,  VEGIE,  and  VCTM  predict  ter¬ 
rain  background  radiance.  The  target  and  background  radiance 
contrast  depends  on  the  relative  values  of  the  corresponding 
surface  emissivities ,  as  well  as  the  surface  temperature  val¬ 
ues.  Therefore,  the  CEES  allows  the  user  to  select  the  value 
of  the  target  surface  emissivity  that  minimizes  the  target- 
background  thermal  contrast  and  hence  minimizes  the  target 
detection  probability  for  any  time  of  day  and  ambient  weather 
conditions. 


Recommendations 


197.  This  study  of  IR,  radar,  and  laser  target  detection  models  and  the 
associated  battlefield  aids  has  produced  the  following  recommendations: 

&.  An  experimental  study  should  be  undertaken  to  evaluate  the 
terrain  signature  models  and  battlefield  aids  described  in 
this  report . 

b .  Signature  models  for  high-value  targets  should  be  improved  and 
embellished. 

c.  Methods  of  improving  target  detection  under  severe  weather 
conditions  should  be  studied. 
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